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ABSTRACT 

We measure the star formation efRciency (SFE), the star formation rate per unit gas, in 23 nearby 
galaxies and compare it to expectations from proposed star formation laws and thresholds. We use 
H I maps from THINGS and derive H2 maps from CO measured by HERACLES and BIMA SONG. 
We estimate the star formation rate by combining GALEX FUV maps and SINGS 24/^m maps, 
infer stellar surface density profiles from SINGS 3.6/im data, and use kinematics from THINGS. We 
measure the SFE as a function of: the free-fall and orbital timescales; midplane gas pressure; stability 
of the gas disk to collapse (including the effects of stars); the ability of perturbations to grow despite 
shear; and the ability of a cold phase to form. In spirals, the SFE of H2 alone is nearly constant at 
5.25 ± 2.5 X 10""'^° yr~^ (equivalent to an H2 depletion time of 1.9 x 10^ yr) as a function of all of 
these variables at our 800 pc resolution. Where the ISM is mostly H I, on the other hand, the SFE 
decreases with increasing radius in both spiral and dwarf galaxies, a decline reasonably described by an 
exponential with scale length 0.2-0.25 r25- We interpret this decline as a strong dependence of CMC 
formation on environment. The ratio of molecular to atomic gas appears to be a smooth function 
of radius, stellar surface density, and pressure spanning from the H2-dominated to H I-dominated 
ISM. The radial decline in SFE is too steep to be reproduced only by increases in the free-fall time 
or orbital time. Thresholds for large-scale instability suggest that our disks are stable or marginally 
stable and do not show a clear link to the declining SFE. We suggest that ISM physics below the 
scales that we observe — phase balance in the H I, H2 formation and destruction, and stellar feedback 
— governs the formation of GMCs from H I. 

Subject headings: galaxies: evolution — galaxies: ISM — radio lines: galaxies — stars: formation 



1. INTRODUCTION 

In nearby galaxies, the star formation rate (SFR) is 
observed to correlate spatially with the distribution of 
neutral gas, at least to first order. This is observed using 
a variety of SFR and gas tracers, but the quantitative re- 
lationship between the two remains poorly understood. 
Although it is common to relate SFR to gas surface den- 
sity via a power law, the relationship is often more com- 
plex. The same surface density of gas can correspond 
to dramatically different SFRs depending on whether it 
is found in a spiral or irregular galaxy or in the inner 
or outer part of a galactic disk. Such variations have 
spurred suggestions that the local potential well, pres- 
sure, coriolis forces, chemical enrichment, or shear may 
regulate the formation of stars from the neutral interstel- 
lar medium (ISM). 

In this paper, we compare a suite of proposed star for- 
mation laws and thresholds to observations. In this way, 
we seek to improve observational constraints on theo- 
ries of galactic-scale star formation. Such theories are 
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relevant to galaxy evolution at all redshifts, but must 
be tested mainly in nearby galaxies, where observations 
have the spatial resolution and sensitivity to map star 
formation to local conditions. An equally important goal 
is to calibrate and test empirical star formation recipes. 
In lieu of a strict theory of star formation, such recipes 
remain indispensable input for galaxy modeling, partic- 
ularly because star formation takes place mostly below 
the resolution of cosmological simulations. This requires 
the implementation of "subgrid" models that map local 
conditions to the SFR (e.g., Springel & Hernquist 2003). 

Our analysis is based on the highest quality data avail- 
able for a significant sample of nearby galaxies: H I maps 
from The H I Nearby Galaxy Survey (THINGS, Walter et 
al. 2008), far ultraviolet (FUV) maps from the GALEX 
Nearby Galaxies Survey (Gil de Paz et al. 2007), infrared 
(IR) data from the Spitzer Infrared Nearby Galaxies Sur- 
vey (SINGS, Kennicutt et al. 2003), CO 1 ^ maps from 
the BIMA Survey of Nearby Galaxies (BIMA SONG, 
Heifer et al. 2003) and CO 2 ^ 1 maps from the HERA 
CO-Line Extragalactic Survey (HERACLES Leroy et al. 
2008). This combination yields sensitive, spatially re- 
solved measurements of kinematics, gas surface density, 
stellar surface density, and SFR surface density across 
the entire optical disks of 23 spiral and irregular galax- 
ies. 

The topic of star formation in galaxies is closely linked 
to that of giant molecular cloud (CMC) formation. In 
the Milky Way, most star formation takes place in 
GMCs, which are predominantly molecular, gravitation- 
ally bound clouds with typical masses ~ 10^ - 10^ M© 
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(Blitz 1993). Similar clouds dominate the molecular ISM 
in Local Group galaxies (e.g., Fukui et al. 1999; Engar- 
giola et al. 2003). If the same is true in other galaxies, 
then a close association between GMCs and star forma- 
tion would be expected to be a general feature of our 
data. Bigiel et al. (2008) study the relationship between 
atomic hydrogen (H I), molecular gas (H2); and the SFR 
in the same data used here. Working at a resolution of 
750 pc, they do not resolve individual GMCs, but do 
find that a single power law with index n= 1.0 ± 0.2 re- 
lates H2 and SFR surface density over the optical disks 
of spirals. This suggests that as in the Milky Way, a key 
prerequisite to forming stars is the formation of GMCs 
(or at least H2). 

Bigiel et al. (2008) find no similar trend relating H I 
and SFR. Instead the ratios of Ha-to-H I and SFR-to-H I 
vary strongly within and among galaxies. CMC forma- 
tion, therefore, appears to be a function of local condi- 
tions. Here we investigate this dependence. We focus 
on where the ISM can form gravitationally bound, pre- 
dominantly molecular structures, i.e., the "star forma- 
tion threshold," and investigate how the molecular frac- 
tion of the ISM varies with local conditions. In equilib- 
rium, the fraction of the ISM in GMCs may be set by the 
timescale over which these structures form. Therefore 
we also consider suggested timescales for the formation 
of GMCs and compare them to observations. 

Maps with good spatial coverage and sensitivity are 
critical to distinguish between the various proposed 
thresholds and timescales. Perhaps the key observation 
to test theories of galactic-scale star formation is that 
the star formation per unit gas mass decreases in the 
outer disks of spiral and irregular galaxies (e.g., Ken- 
nicutt 1989; Martin & Kennicutt 2001; Thornley et al. 
2006). The details of this decrease vary with the specifics 
of the observations. For example, Martin & Kennicutt 
(2001) observed a sharp drop in the distribution of H II 
regions, while UV maps suggest a steady decline (Boissier 
et al. 2007), but it is without dispute that the SFR per 
unit gas mass does indeed decline (see also Wong & Blitz 
2002). Maps with good spatial extent contain both re- 
gions where CMC formation proceeds efhciently and re- 
gions where it is suppressed. Including both H I-rich 
dwarf galaxies and Ha-dominated spirals offers a simi- 
lar contrast. 

In §2, we present a set of star formation laws and 
thresholds that we will compare to observations. We 
phrase these in terms of the star formation per unit neu- 
tral gas, which wc call the "star formation efhciency" 
This quantity, the inverse of the gas depletion time, re- 
moves the basic scaling between stars and gas and mea- 
sures how effectively each parcel of the ISM forms stars. 

In §3, wc briefly describe our sample, data, and 
methodology. In order to focus the main part of the 
paper on analysis, we defer most detailed discussion of 
data and methodology to the appendices. 

In §4 we look at how the star formation efficiency 
relates to other basic quantities (§4.1), proposed laws 
(§4.2), and thresholds (§4.3) described in §2. In §5 wc 
analyze and interpret these results. In §6, we illustrate 
our conclusions by comparing predictions for the star for- 
mation efficiency to observations. In §7, we summarize 
our results. 

Appendices A - D contain all the information required 



to reproduce our calculations, including descriptions of 
the data and how we convert from observables to physical 
quantities. We present our data as an electronic table of 
radial profiles described in Appendix E and as maps and 
plotted profiles for each galaxy in Appendix F. 

2. BACKGROUND 

Following, e.g., Kennicutt (1989), we break the topic 
of star formation in galaxies into two parts. Where star 

formation is widespread, we refer to the quantitative re- 
lationship between neutral gas and the SFR as the star 
formation law. To predict the SFR over an entire galac- 
tic disk, it is also necessary to know which gas is actively 
forming stars. This topic is often phrased as the star 
formation threshold, but may be more generally thought 
of as the problem of where a cold phase (n ~ 4-80 cm~^, 
T ^ 50-200 K) or gravitationally bound clouds can form; 
both are thought to be prerequisites to star formation. 
We give a brief background on both laws and thresholds, 
first noting that neither term is strictly accurate: "laws" 
here refer to observed (or predicted) correlations and the 
"threshold" is probably a smooth variation from non-star 
forming to actively star forming gas. 

We cast this discussion in terms of the star forma- 
tion efficiency (SEE). There are many definitions for the 
SEE, but throughout this paper we use the term only 
to refer to the star formation rate surface density per 
unit neutral gas surface density along a line of sight, 
i.e., SFE=SsFR/Sgas with units of yr~^. We will also 
discuss SEE (H2) which refers to the SFR per unit H2 
(Esfr/Sh2), and SEE (H I)^ (Ssfr/Shi). The SEE is 
the inverse of the gas depletion time, the time required 
for present day star formation to consume the gas reser- 
voir. It represents a combination of the real timescale 
for neutral gas to form stars and the fraction of gas that 
ends up in stars, e.g., if 1% of the gas is converted to 
stars every 10^ yr, the SEE= 10^^ yr^^- Because it is 
normalized by Sgas, the SEE is more useful than Ssfr 
alone to identify where conditions are conducive to star 
formation (i.e., where gas is "good at forming stars"). 

As we describe proposed laws (§2.1) and thresholds 
(§2.2), we present quantitative forms for each that can 
be compared to the observed SEE. Table 1 collects these 
expressions, which we compare to observations in §4. 

2.1. Star Formation Laws 

A star formation law should predict the SEE from local 
conditions. Here we describe three proposals for the lim- 
iting timescale over which gas forms stars: the free-fall 
timescale in the gas disk, the orbital timescale, and the 
characteristic timescale for cloud-cloud collisions. We 
also describe proposals that GMCs form stars with a 
fixed SEE and that the midplane gas pressure regulates 
the fraction of the ISM in the molecular phase. We 
present each proposal as a prediction for the SFE in 
terms of observables. These appear together in the upper 
part of Table 1. We expect a successful star formation 
law to reproduce the observed SFE (in practice, com- 
bined with an empirical calibration). 

2.1.1. Disk Free-Fall Time With Fixed Scale Height 

The most common formulation of the star formation 
law is a power law relating gas and star formation (sur- 
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TABLE 1 






Star Formation Laws and Thresholds 




Theory 


Form 


Observablcs 


Star Foriiialioii Laws 



. . . fixed scale height 
. . . variable scale height 

orbital timcscalc 
cloud-cloud collisions 
fixed GMC efiiciency 

pressure and ISM phase 



SFE oc EO^I 



SFE or il^ol oc ^ (l + ^^) 

SFE or R„ 



nol Tnrh ~ 



"("•ral) 



orb 27rrgal 

SFE cx Qias (1 - 0.7 /3) 
SFE = SFE (H2) Tr"^ 

-^mol ^ ( ^gas ( Sgas + ^ E* 



1.2 



^gas 5 ^"3 5 C"* 

''('•gal) 

f(rgal) 

Egas 5 (Tg, (J* 



Star Formation Thresholds 



gravitational instability 

... in the gas disk Qgas = ^ c's . 

... in a disk of gas and stars Qgti — f 2 

competition with shear E 
cold gas phase 



< 1 



ars+gas — \ /~i , i i „ 
^ 2.5 A rr 



Egal > 6.i"M0 pC-2 /0-3 Z-0-3 70-23 



Egas, 0-g, -"(j-gal) 

Egas, S,, o-g, cr,, D(rgal) 

Egas, O-g, ^(rgal) 
Egas, E*, Zy I 



face) densities following Schmidt (1959, 1963). Kenni- 
cutt (1989, 1998a) calibrated this law in its observable 
(surface density) form. Averaging over the star-forming 
disks of spiral and starburst galaxies, he found 



^SFR 



CX S 



1.4 



(1) 



often referred to as the "Kcnnicutt Schmidt law." 

The exponent in Equation 1, n « 1.5, can be ap- 
proximately explained by arguing that stars form with 

a characteristic timcscale equal to the free fall time in 
the gas disk, which in turn depends inversely on the 
square root of the gas volume density, Tg oc Pg^^^ (c-g-, 
Madore 1977). For a fixed scale height pg^s « Egas and 
SsFR oc Sg^g. The first star formation law that we con- 
sider is thus 



SFE (X Eli , 
which is approximately Equation 1. 



(2) 



2.1.2. Disk Free-Fall Time With Variable Scale Height 

If the scale height is not fixed, but instead set by hy- 
drostatic equilibrium in the disk, then 



1 



Tft CX 



(3) 



where o-g and it*,^ are the (vertical) velocity dispersions 
of gas and stars, Egas and are the surface densities of 
the same, and Pmp.gas is the midplane gas density. Equa- 
tion 3 combines the expression for midplane density from 
Krumholz & McKee (2005, their Equation 34) and mid- 
plane gas pressure from Elmegreen (1989, his Equation 
11, used to calculate (j)p). The second star formation law 
that we consider is 



SFEoc 



1 



0.5 



(4) 



which incorporates variations in the scale height and thus 
gas volume density with a changing potential well. 



2.1.3. Orbital Timescale 

It is also common to equate the timcscalc for star 
formation and the orbital timescale (e.g.. Silk 1997; 
Elmegreen 1997). Kcnnicutt (1998a) and Wong & Blitz 
(2002) found that such a formulation performs as well as 
Equation 1. In this case 



SFEocr- = i^ = ^. 

2n 27rrgai 

where f (rgai) is the rotational velocity at a galactocentric 
radius rgai and is the corresponding angular velocity. 

2.1.4. Cloud-Cloud Collisions 

Tan (2000) suggested that the rate of collisions between 
gravitationally bound clouds sets the timescale for star 
formation so that 

SFE CX T-l Qg-i (1 - 0.7^) . (6) 

where Qgas, defined below, measures gravitational insta- 
bility in the disk and /3 = dlog i;(rgai)/rflog rgai is the 
logarithmic derivative of the rotation curve. The depen- 
dence on f3 reflects the importance of galactic shear in 
setting the frequency of cloud-cloud collisions. In the 
limit /3 = (a flat rotation curve) this prescription re- 
duces to essentially Equation 5; for /? = 1 (solid body 
rotation) the SFE is depressed by the absence of shear. 

2.1.5. Fixed GMC Efficiency 

If the SFE of an individual GMC depends on its intrin- 
sic properties and if these properties are not themselves 
strong functions of environment or cloud formation, then 
we expect a fixed SFR per unit molecular gas, SFE (H2). 
Krumholz & McKee (2005) posited such a case, arguing 
that the SFE of a GMC depends on the free-fall time 
in the cloud, itself only a weak function of cloud mass in 
the Milky Way (Solomon et al. 1987). Bigiel et al. (2008) 
found support for this idea. Studying the same data used 
here, they derived a linear relationship between Eh2 and 
^SFR on scales of 750 pc. 

SFE (II2) is likely to appear constant if: the scaling re- 
lations and mass spectrum (i.e., the intrinsic properties) 



(5) 
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of GMCs arc approximately universal, the gas pressure 
is low enough that GMCs arc largely decoupled from the 
rest of the ISM, individual resolution elements contain 
at least a few GMCs, and the properties of a cloud regu- 
late its ability to form stars (§5.1 and Bigiel et al. 2008). 
This is the fifth star formation law that we consider, that 
star formation in spiral galaxies occurs mostly in GMCs 
and that once such clouds are formed, they have approx- 
imately uniform properties so that 

SFE (H2) = constant , (7) 

which we can convert to the SFE of the total gas given 
^^moi = Sh2/Shi, the ratio of H2 to H I gas. Then 



SFE = SFE (H2) 



Rn 



(8) 



-Rmol + 1 

or if we measure only Shi (as is the case in dwarfs) , then 
SFE (HI) = SFE (H2) R^oi. 

The balance between CMC /H2 formation and destruc- 
tion will set Rmoi = Sh2/Shi. If GMCs with fixed 
lifetime form over a free fall time or orbital time then 
-Rmoi oc r^^ or -Rmoi oc T~l^ (§5.4), which we have noted 
in Tabic 1. Combined with Equation 8, an expression for 
-Rmoi predicts the SFE. 

2.1.6. Pressure and Phase of the ISM 

Wong & Blitz (2002), Blitz & Rosolowsky (2004), and 
Blitz & Rosolowsky (2006) explicitly consider i?inoi- Fol- 
lowing Elmegreen (1989) and Elmcgrcen & Parravano 
(1994), they identify pressure as the critical quantity 
that sets the ability of the ISM to form H2 . They show 
that the midplane hydrostatic gas pressure, Ph, corre- 
lates with this ratio in the inner parts of spiral galaxies. 

Pressure, which is directly proportional to the gas 
volume density, should affect both the rate of H2 for- 
mation/destruction and the likelihood of a gravitation- 
ally unstable overdensity condensing out of a turbulent 
ISM (Elmegreen 1989; Elmegreen & Parravano 1994). 
Elmegreen (1989) gives the following expression for Ph, 



Ph 



— U A,gas 



(7* 



(9) 



and Elmegreen (1993) predicted that the fraction of gas 
in the molecular phase depends on both Ph and the in- 
terstellar radiation field, j, via Pmoi oc P^'^ If 
5^SFR OC Sh2 and wc make the simple assumption that 
j oc SsFR then Elmegreen (1993) predicts 



Pmoi « Ph or Sh2 = ^mPh , (10) 
which combines with Equation 8 to predict the SFE. 

Wong & Blitz (2002) and Blitz & Rosolowsky (2006) 
found observational support for Equation 10. Using a 
modified Equation 9 appropriate where > Sgas, Blitz 
& Rosolowsky (2006) fit a power law of the form 



R 



mol 



Sh2 

Shi 



(11) 



finding Pq = 4.3 x lO** cm K, the observed pressure 
where the ISM is equal parts H I and H2, and a best- 
fit exponent a = 0.92. Wong & Blitz (2002) found 
a — 0.8. Robertson & Kravtsov (2008) recently found 
support from simulations for a 0.9. 



2.2. Star Formation Thresholds 

We have described suggestions for the efficiency with 

which gas form stars, but not whether gas forms stars. 
A "star formation threshold" is often invoked to accom- 
pany a star formation law. This is a criterion designed 
to address the question "which gas is actively forming 
stars?" or "where can the ISM form gravitationally 
bound, molecular clouds?" and proposed thresholds have 
mostly focused on the existence of gravitational or ther- 
mal instability in the gas disk. 

A common way to treat the issue of thresholds is to 
formulate a critical gas surface density, Scrit, that is a 
function of local conditions — kinematics, stellar surface 
density, or metallicity. If Sgas is below Ecrit, star for- 
mation is expected to be suppressed; we refer to such 
regions as "subcritical." Where the gas surface density 
is above the critical surface density, star formation is ex- 
pected to be widespread. We refer to such regions as 
"supercritical." 

In practice, wc expect to observe a drop in the SFE 
associated with the transition from super- to subcritical. 
We do not necessarily expect SFE = in subcritical 
regions. Even with excellent resolution, a line of sight 
through a galaxy probes a range of physical conditions. 
At our working resolution of 400 800 pc, each resolution 
element encompasses a wide range of local conditions. 
Within a subcritical resolution element, star formation 
may still occur in isolated pockets that locally meet the 
threshold criterion. 

Expressions for star formation thresholds are collected 
in the lower part of Table 1. 

2.2.1. Gramtational Instability 

Kcnnicutt (1989), Kcnnicutt (1998a), and Martin & 
Kennicutt (2001) argued that star formation is only 
widespread where the gas disk is unstable against large 
scale collapse. Following Toomre (1964), the condition 
for instability in a thin gas disk is 



Q 



gas 



TT G S„ 



< 1 . 



(12) 



where ag is the gas velocity dispersion, G is the gravita- 
tional constant, and k is the epicyclic frequency, calcu- 
lated via 



1.41 



^('■gal) 



''gal 



(13) 



where (3 = dXog w(rgai)/dlog rgai. 

Martin & Kcnnicutt (2001) found that H II regions are 
common where Sgas exceeds a critical surface density de- 
rived following Equation 12, 

^crit.Q = OlQ 



G 



(14) 



In regions where Sgas is above this threshold, gas is un- 
stable against large scale collapse, which leads to star 
formation. Below the threshold, Coriolis forces counter- 
act the self-gravity of the gas and suppress cloud/star 
formation. The factor aq is an empirical calibration, 
the observed average value of l/Qgas at the star for- 
mation threshold. For an ideal thin gas disk, the con- 
dition for gas to be unstable to collapse is ag > 1. 
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At the edge of star forming disks, Kennicutt (1989) 
found QQ =- 0.63 and Martin & Kennicutt (2001) found 
aq = 0.69 (Qgas ~ 1.5). 

Kennicutt (1989) and Martin & Kennicutt (2001) men- 
tion the influence of stars as a possible cause for Qgas > 1 
at the star formation threshold. Hunter et al. (1998a) 
present an in-depth discussion of how several factors in- 
fluence ag, e.g., stars and viscosity lower it, while the 
thickness of the gas disk raises it. Kim & Ostriker 
(2001, 2007) argue based on simulations that the ob- 
served threshold corresponds to the onset of nonlinear, 
non-axisymmetric instabilities. Schayc (2004) and dc 
Blok & Walter (2006) suggest a different explanation, 
that OLQ ^ \ because Og has been systematically mis- 
handled; they point out that Og measured from 21-cm 
emission will overestimate the true velocity dispersion of 
gas in a cold phase. 

The stellar potential well may substantially affect the 
stability of the gas disk. Rafikov (2001) extended work by 
Jog & Solomon (1984) to provide a straightforward way 
to calculate the instability of a gas disk in the presence 
of a coUisionless stellar disk. Rafikov defined 



Qstars — 



(7 ^ f tx 



(15) 



where (t,^^ is the (radial) velocity dispersion of stars and 
is the stellar mass surface density. The condition for 
instability in the gas disk is 



^stars-|-ga 



Qstars 1 + 9^ Qg 



R 



1 + q'E? 



> 1 



where q = fc(T,.r/K, with k the wavenumbcr of the insta- 
bility being considered, and R = (7g/<7*,r- The minimum 
value of Qstars-i-gas indicates whether the gas disk is un- 
stable to large scale collapse. In our sample, typical val- 
ues of q correspond to wavelengths A = 27r/fc « 1-5 kpc. 

Hunter ct al. (1998a) and BHtz & Rosolowsky (2004) 
observed strong correlations between star and GMC for- 
mation and the distribution of stars, consistent with stel- 
lar gravity playing a key role in star formation. Yang et 
al. (2007) recently showed that Qstars+gas does an ex- 
cellent job of predicting the location of star formation 
in the Large Magellanic Cloud and Boissier et al. (2003) 
showed that including stars improves the correspondence 
between Q and star formation in disk galaxies. Li et al. 
(2005, 2006) found the same results from numerical sim- 
ulations of disk galaxies, i.e., that stability against large 
scale collapse depends critically on the stellar potential 
well, with star formation where Qstars-Fgas ^ 1-6. 

2.2.2. Galactic Shear 

Motivated by the failure of the Toomre Qgas threshold 
in dwarf irregular galaxies, Hunter et al. (1998a) sug- 
gested that collecting the material for cloud formation 
may be easier than implied by Qgas, e-g-, through the 
aid of magnetic fields (sec also Kim & Ostriker 2001). 
They hypothesize that the destructive influence of galac- 
tic shear may instead limit where GMCs can form and 
describe a threshold that depends on the ability of clouds 
to form in the time allowed by shear. 

This threshold is based on the local shear rate, de- 
scribed by Oort's A constant 



A = -0.5 Tgal 

Substituting = w(rgai)/rgai. 



dr. 



gal 



(17) 



A = 0.5 r^-^"/^-')Uo.5^ (1-/3) 

(18) 



Then the threshold has the form 



-'crit.A 



aA cTg A 

TT G 



(19) 



Hunter et al. (1998a) suggest aA = 2.5, but this nor- 
malization for Scrit.A is relatively uncertain. The value 
chosen by Hunter et al. (1998a) corresponds to pertur- 
bations growing by a factor of 100 during the time 
allowed by shear, which roughly matches both the sur- 
face density contrast between Shi and a GMC and the 
condition Qgas S 1 where (iw(rgai)/(irgai = 0. 

The practical advantage of shear over Qgas is that shear 
is low in dwarf galaxies and the inner disks of spiral galax- 
ies {(3 = 1 for solid body rotation), both locales where 
widespread star formation is observed. In the outer disks 
of spiral galaxies — where star formation cutoffs are ob- 
served — rotation curves tend to be flat {f3 = 0) so that 
^critjA and Ecrit.Q rcducc to the same form. 

2.2.3. Formation of a Cold Phase 

The very long time needed to assemble a massive GMC 
from coagulation of smaller clouds suggests that most 
GMCs in galaxy disks form "top down" (e.g., McKee & 
Ostriker 2007). However this does not necessarily require 
that the whole gas disk to be unstable. Where cold H I is 
abundant, the lower velocity dispersion associated with 
this phase may render the ISM locally unstable (Schaye 
2004), leading to the formation of GMCs and stars. 

Therefore, instead of large-scale gravitational instabil- 
ity or cloud destruction by shear, the ability to form a 
cold neutral medium (McKee & Ostriker 1977; Wolflre et 
al. 2003) may regulate GMC formation. Schayc (2004) 
argues based on modeling that near the cutoffs observed 
by Martin & Kennicutt (2001) gas becomes mostly cold 
H I and H2, (Tg drops accordingly, and Q becomes < 1 
in the cold gas. In a similar vein, Elmegreen & Par- 
ravano (1994) suggest that the star formation efliciency 
in the outer parts of galaxies drops because the pres- 
sure becomes too low to allow a cold phase to form even 
given perturbations, e.g., from supernova shocks. Braun 
(1997) found support for this idea using 21-cm observa- 
tions; he associated networks of high surface brightness 
filaments with cold H I and showed that these filaments 
are pervasive across the star forming disk, but become 
less common at large radii (though work on THINGS by 
Usero et al. 2008, calls this result into question). 

Schaye (2004) modeled the ISM to estimate where the 
average temperature drops to « 500 K, the molecular 
fraction reaches w 10~^, and Qgas ~ 1; good indicators 
that cold H I is common and H2 formation is efficient. 
These all occur where Sgas exceeds 



6.1 



-'S04 



Jf 



0.3 



M© pc-2 •'s \0.1Z(, 



-0.3 



0.23 



106 cm-2 s-1 J 
(20) 



6 



Leroy et al. 



where /g w Egas/(Sgas + S,) is the fraction of mass in gas 
(we assume a two-component disk) , Z is the mctallicity of 
the ISM, and / is the flux of ionizing photons. Sso4 also 
depends on the ratio of thermal to turbulent pressure and 
higher order terms not shown here. Schaye (2004) selects 
fiducial values to match those expected in outer galaxy 
disks, but concludes that the influence of Z, /g, and the 
radiation field is relatively small. Most reasonable values 
yield Sso4 « 3 - 10 M© pc'^. 

Schaye (2004) argues that a simple column density 
threshold may work as well as dynamical thresholds. 
This agrees with the observation by, e.g., Skillman (1987) 
and de Blok & Walter (2006) that a simple H I column 
density threshold does a good job of predicting the loca- 
tion of star formation in dwarf irregulars. This thresh- 
old, Ehi ~ 10 M© pc^^, also corresponds to the surface 
density above which H I is observed to saturate (Mar- 
tin & Kennicutt 2001; Wong & Blitz 2002; Bigiel et al. 
2008); that is, gas in excess of this surface density in 
spiral galaxies is in the molecular phase. 

3. DATA 

The right hand column of Table 1 lists the observables 
required to evaluate each law or threshold. We require 
estimates of: the surface density of atomic gas (Shi), 
molecular gas (Sh2), star formation rate (Ssfr), and 
stellar mass (E*), the velocity dispersions of gas and stars 
(cgas and (7*), and the rotation curve (v(r^^\)). Estimates 
of the metallicity await future work. 

3.1. The Sample 

We assemble maps and radial profiles of the necessary 
quantities in 23 nearby, star-forming galaxies that we 
list in order of increasing stellar mass in Table 2. These 
are galaxies for which we could compile the necessary 
data, which means the overlap of THINGS, SINGS, the 
GALEX NGS, and (for spirals) either BIMA SONG or 
HERACLES. 

We work with two subsamples: 11 H I-dominated, low- 
mass galaxies and 12 large spiral galaxies. In Table 2, the 
galaxies that we classify "dwarf galaxies" lie above the 
horizontal dividing line. These have rotation velocities 
Vrot ^ 125 km s"^, stellar masses < 10^° M©, and 
Mb > —20 mag. The galaxies that we label "spirals'' 
lie below the dividing line and have v^ot ^ 125 km s~^, 
M* > 10^° M©, and Mb < -20 mag. 

This division allows us to explore two distinct regimes 
in parallel. Compared to their larger cousins, dwarf 
galax;ies have low metallicities, intense radiation fields, 
lower galactic shear, and weak or absent spiral structure. 
Metallicity, in particular, should have a strong effect on 
the thermal balance of the ISM. In lieu of direct mea- 
surements, separating the sample in this way allows us 
to assess its impact. 

We treat the two subsamples slightly differently in two 
ways. First, we place data for spirals at a common spa- 
tial resolution of 800 pc and data for dwarf galaxies at 
400 pc. The spirals in our sample are farther away than 
the dwarf galaxies with larger physical radii, and this 
approach ensures a good number of resolution elements 
across each galaxy and a fairly uniform angular resolu- 
tion of ~ 20" (see Table 2). 

Second, we use CO maps combined with a constant 
CO-to-H2 conversion factor, Xco, to derive Eh2 in spi- 



TABLE 2 
Sample Galaxies 



Galaxy^ 


TD b 

Kes. 


GO 


Rotation 


Also in 




(") 




Curve"^ 


sample of* 




19 




dB 




Ho I 


21 




T 




Ho H 


24 




T 




IC 2574 


21 




dB 




NGC 4214" 


28 




T 




NOC 2Q7fi 


23 




dB 




NGC 4449° 


20 








NGC 3077'' 


22 








NOC 77qs 


21 




dB 




NGC 925 


9 




dB 


12 4 


NGC 2403 


26 




dB 


12 4 


NGC 628 


23 


HERACLES 


T 


1, 2 


NGC 3198 


12 


HERACLES 


dB 




NGC 3184 


15 


HERACLES 


T 




NGC 4736 


35 


HERACLES 


dB 


1, 2, 3, 5 


NGC 3351 


16 


HERACLES 


T 




NGC 6946 


28 


HERACLES 


dB 


2 


NGC 3627 


18 


BIMA SONG 


dB 


5 


NGC 5194 


21 


BIMA SONG 


T 


2, 4, 5 


NGC 3521 


15 


HERACLES 


dB 


5 


NGC 2841 


12 


HERACLES 


dB 


1, 2 


NGC 5055 


16 


HERACLES 


dB 


2. 3. 5 


NGC V:5:-!l 


11 


ilEl^ACLES 


dii 


2. r, 



^'In order of increasing stellar mass. 

''Angular resolution to match working spatial resolution 
in the subsample, 400 pc for dwarf galaxies and and 800 pc 
for spirals. 

Rotation curve data: dB = de Blok et al. (2008); T = 
only THINGS first moment (Walter et al. 2008) 

dl: Kennicutt (1989); 2: Martin & Kennicutt (2001); 3: 
Wong & Blitz (2002); 4: Boissier et al. (2003); 5: Blitz & 
Rosolowsky (2006) 

•^IR data from Spitzer archive (not SINGS). 

rals, while we treat the molecular gas content of dwarf 
galaxies as unknown (see Appendix A. 3). CO emission 
in very low mass galax;ies is usually weak or not detected 
(e.g., Taylor et al. 1998; Leroy et al. 2005, and see Table 
4) and its interpretation is confused by potential varia- 
tions in Xqo ■ Because dwarf galaxies lack H2-filled H I 
depressions like those observed in the centers of spirals, 
we expect Ehi to at least capture the basic morphology of 
the total gas. Although we do not measure Sh2 in dwarf 
galaxies, we consider our results in light of the possibility 
of an unseen reservoir of molecular gas (§5.3). 

3.2. Data to Physical Quantities 

Appendices A D explain in detail how we translate 
observables into physical quantities. Here and in Table 
3 we summarize this mapping. 

Atomic Hydrogen Surface Density (Appendix A): We 
derive atomic gas mass surface density, Ehi, from 21- 
cm line integrated intensity maps obtained by Walter 
et al. (2008) as part of the THINGS survey using the 
Very Large Array ^. Ehi if' corrected for inclination and 
includes a factor of 1.36 to account for helium. 

Molecular Hydrogen Surface Density (Appendix A): 
In spirals, we estimate the molecular gas mass surface 
density, Sh2, from CO line emission. For 10 galaxies 

^ The VLA is operated by the National Radio Astronomy Ob- 
servatory, which is a facility of the National Science Foundation 
operated under cooperative agreement by Associated Universities, 
Inc. 
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TABLE 3 
Data to Physical Quantities 



Quantity 



Observation 



Survey 



Reference 



Key Assumptions 



Shi 

Sh2 (spirals only) 

Unobscured SgpR 
Embedded Ssfr 

IviiieiiiaLics 



21-cm line 
CO 2 -» 1 
CO 1 -»0 
FUV 
24/im 
3. 6/ /111 
21-ciii liiie 



THINGS 

HERACLES 

BIMA SONG 

GALEX NGS 

SINGS 

SINGS 

TlllNCiS 



Walter et al. (2008) 
Leroy et al. (2008) 
Heifer et al. (2003) 
Gil de Paz et al. (2007) 
Kennicutt et al. (2003) 
Krimirntt ct al. (2003) 
de Bluk el al. (2008) 



fixed line ratio, CO-to-H2 conversion 
fixed CO-to-H2 conversion 



Tf = 0..5 M-: /L -. 

tihiipk} iuiicLioiiai liL; lixcxl (Tga 



we use data from HERACLES, a large program at the 
IRAM^ 30-m telescope (Leroy et al. 2008) that used the 
HERA focal plane array (Schuster et al. 2004) to map a 
subsample of THINGS' in the CO J = 2 ^ 1 line. For 
NGC 3627 and NGC 5194, we use J = 1 line maps 
from the BIMA SONG survey (Heifer et al. 2003). 

We convert from CO line intensity to Sh2 as- 
suming a constant C0-to-H2 conversion factor ap- 
propriate for the solar neighborhood, Xqo = 2 x 
10^° cm~^ (K km s~^)~^, and a fixed line ratio /co(2 
1) = 0.8 /co(l ~* 0), typical of the disks of spiral galax- 
ies. We correct for the effects of inclination and include 
a factor of 1.36 to reflect the presence of helium . 

Galactic Rotation (Appendix B): We fit a simple func- 
tional form to the high quality rotation curves derived 
from THINGS by de Blok et al. (2008) and the THINGS 
first moment maps (Walter et al. 2008). These fits yield 
smooth, well-behaved (analytic) derivatives and match 
the observations well. Two galaxies (NGC 3077 and 
NGC 4449) have complex velocity fields that require sub- 
stantial effort to interpret and we omit them from anal- 
yses requiring kinematics. 

Gas Velocity Dispersion (Appendix B): We assume a 
fixed gas velocity dispersion, dgas = 11 km s""'^, a value 
motivated by the THINGS second moment maps. 

Stellar Velocity Dispersion (Appendix B): We estimate 
the vertical stellar velocity dispersion, ct*,z, from hydro- 
static equilibrium, the assumption of an isothermal disk, 
and an estimated (radially invariant) stellar scale height. 
We derive this scale height for each galaxy from our mea- 
sured stellar scale length and an average flattening ratio 
for disk galaxies. We take the vertical and radial velocity 
dispersions to be related by (t,^z = 0.6 u*,,-- 

Stellar Surface Density (Appendix C): We estimate 
the stellar surface density, S,, from Spitzer 3.6/Ltm maps, 
mostly from SINGS (Kennicutt et al. 2003). To avoid 
contamination by hot dust and foreground stars, we con- 
struct radial profiles only, using the median 3.6/xm inten- 
sity in each tilted ring. We convert from 3.6/im intensity 
to via an empirical K-to-',i.6fj,m calibration and adopt 
a fixed K-h&nd mass-to-light ratio, = 0.5 Mq/Lq^k- 

Star Formation Rate Surface Density (Appendix D): 
We combine FUV and 24/zm maps to derive maps of 
EgpR,; giving us a tracer sensitive to both exposed and 
dust-embedded star formation. The FUV data come 
from the GALEX Nearby Galaxies Survey (Gil de Paz 
et al. 2007) and the 24/im maps are part of SINGS. Be- 
cause this precise combination of data is new, Appendix 
D includes an extended motivation for how we convert 
intensity to SgpR- 

8 IRAM is supported by CNRS/INSU (France), the MPG (Ger- 
many) and the IGN (Spain) 



3.3. Properties of the Sample 

Table 4 compiles the integrated properties of each 
galaxy in our sample. Columns (1) - (7) give basic 
parameters adopted from other sources: the name of 
the galaxy; the distance, inclination, and position an- 
gle (Walter et al. 2008, except that we adopt i = 20° in 
M 51); and the morphology, B-band isophotal radius at 
25 mag arcsec"^ (''25), and B-band absolute magnitude 
from LEDA (Prugniel & Heraudeau 1998). Columns (8) 
and (9) give Uflat and Zfiat, the free parameters for our 
rotation curve fit (Appendix B); from these two param- 
eters one can calculate v (rgai) and /?. Columns (10) - 
(13) give the total stellar mass, H 1 mass, H2 mass and 
SFR from integrating our data within 1.5 r25. 

Columns (14) (17) give scale lengths derived from 
exponential fits to the S*, SgpR, and Sh2 (CO) radial 
profiles. The stellar scale lengths match those found by 
Tamburro et al. (2008) with 15% scatter; they are ~ 10% 
shorter than those found by Regan et al. (2001), with 
RMS scatter of 20%. Our CO scale lengths are taken 
from Leroy et al. (2008); these are ~ 30% shorter than 
those of Regan et al. (2001) on average. 

3.4. Methodology 

We work with maps of Shi, 5]h2 and EgpR on the 
THINGS astrometric grid. All data are placed at a com- 
mon spatial resolution, 400 pc for dwarf galaxies and 
800 pc for spirals; when necessary, we use a Gaussian ker- 
nel to degrade our data to this resolution. The convolu- 
tion occurs before any deprojcction and may be thought 
of as placing each subsample at a single distance. Radial 
profiles of these maps and appear in Appendix E. 

Using these data, we compute each quantity in Table 1 
for each pixel inside 1.2 r25 and derive radial profiles over 
the same range following the methodology in Appendix 
E. Because we measure and i'(rgai) only in radial pro- 
file, these maps are often a hybrid between radial profiles 
and pixel-by-pixel measurements. 

In §4 - 6, we analyze the combined data set for the 
two subsamples and avoid discussing results for individ- 
ual galaxies. We refer readers interested in individual 
galaxies to the Appendices. Appendix E gives our radial 
profile data and the atlas in Appendix F shows maps of 
5^Hi, 5]h2, total gas, unobscured Ssfr, dust-embedded 
Ssfr , and total Ssfr, as well as profiles of the quanti- 
ties in Table 1. 

In keeping with our emphasis on the combined datasot, 
we default to quoting the mean and la scatter when we 
give uncertainties in parameters derived from the ensem- 
ble of galaxies (we usually estimate the scatter using the 
median absolute deviation to reduce sensitivity to out- 
liers). We prefer this approach to giving the uncertainty 



8 



Leroy et al. 



TABLE 4 

Properties of Sample Galaxies 



Galaxy Dist. 

(Mpc) 


i 

D 


FA 

(°) 


— ; 

Morpn. 




Mb 
(mag) 


r25 

(kpc) 


(km s ^) 


— 

'flat 

(kpc) 


— — 

log Mt 

(M©) 




log Mhi 

(Me) 




log Mh2 

(M©) 


SFR 
(M0 yr-i) 




(kpc) 




«SFR 

(kpc) 




'CO 
(kpc) 




4.3 


66 


230 


Irr 


-14.4 


1.2 


50 


2.0 


7.1 


8.7 


< 6.8 


0.005 


0.8 


1.0 




Ho 1 


3.8 


12 


50 


Irr 


-14.9 


1.8 


53 


0.4 


7.4 


8.3 


< 7.2 


0.009 


0.8 


1.2 




Ho n 


3.4 


41 


177 


Irr 


-16.9 


3.7 


36 


0.6 


8.3 


8.9 


< 7.6 


0.048 


1.2 


1.3 




IC 2574 


4.0 


53 


56 


Irr 


-18.0 


7.5 


134 


12.9 


8.7 


9.3 


< 7.9 


0.070 


2.1 


4.8 




NGG 4214 


2.9 


44 


65 


Irr 


-17.4 


2.9 


57 


0.9 


8.8 


8.7 


~7.0 


0.107 


0.7 


0.5 




NGG 2Q7fi 


3.6 


65 


335 


Sc 


-17.8 


3.8 


92 


1.2 


9.1 


8.3 


7.8 


0.087 


0.9 


0.8 


1.2 


NGG 444Q 


4.2 


60 


230 




-19.1 


2.8 






9.3 


9.2 




0.371 


0.9 


0.8 




NGG Sn77 


3.8 


46 


45 


Sd 


-17.7 


3.0 






9.3 


9.1 




0.086 


0.7 


0.3 




NGG 77qS 


3.9 


50 


290 


Scd 


-18.7 


6.0 


115 


1.5 


9.5 


9.1 




0.235 


1.3 


1.3 




NGC 2403 


3.2 


63 


124 


SBc 


-19.4 


7.3 


134 


1.7 


9.7 


9.5 


7.3 


0.382 


1.6 


2.0 


1.9 


NGC 0925 


9.2 


66 


287 


SBcd 


-20.0 


14.2 


136 


6.5 


9.9 


9.8 


8.4 


0.561 


4.1 


4.1 




NGC 0628 


7.3 


J 


20 


Sc 


-20.0 


10.4 


217 


0.8 


10.1 


9.7 


9.0 


0.807 


2.3 


2.4 


2.4 


NGC 3198 


13.8 


72 


215 


SBc 


-20.7 


13.0 


150 


2.8 


10.1 


10.1 


8.8 


0.931 


3.2 


3.4 


2.7 


NGC 3184 


11.1 


16 


179 


SBc 


-19.9 


11.9 


210 


2.8 


10.3 


9.6 


9.2 


0.901 


2.4 


2.8 


2.9 


NGC 4736 


4.7 


41 


296 


Sab 


-20.0 


5.3 


156 


0.2 


10.3 


8.7 


8.6 


0.481 


1.1 


0.9 


0.8 


NGC 3351 


10.1 


41 


192 


SBb 


-19.7 


10.6 


196 


0.7 


10.4 


9.2 


9.0 


0.940 


2.2 


1.8 


2.5 


NGC 6946 


5.9 


33 


243 


SBc 


-20.9 


9.8 


186 


1.4 


10.5 


9.8 


9.6 


3.239 


2.5 


2.7 


1.9 


NGC 3627 


9.3 


62 


173 


SBb 


-20.8 


13.9 


192 


1.2 


10.6 


9.0 


9.1 


2.217 


2.8 


1.9 


2.2 


NGC 5194 


8.0 


20 


172 


SBc 


-21.1 


9.0 


219 


0.8 


10.6 


9.5 


9.4 


3.125 


2.8 


2.4 


2.3 


NGC 3521 


10.7 


73 


340 


SBbc 


-20.9 


12.9 


227 


1.4 


10.7 


10.0 


9.6 


2.104 


2.9 


3.1 


2.2 


NGC 2841 


14.1 


74 


153 


Sb 


-21.2 


14.2 


302 


0.6 


10.8 


10.1 


8.5 


0.741 


4.0 


5.3 




NGC 5055 


10.1 


59 


102 


Sbc 


-20.6 


17.4 


192 


0.7 


10.8 


10.1 


9.7 


2.123 


3.2 


3.1 


3.1 


NGC 7331 


14.7 


76 


168 


SAb 


-21.7 


19.6 


244 


1.3 


10.9 


10.1 


9.7 


2.987 


3.3 


4.5 


3.1 



^Unless noted log Mh2 comes from HERACLES Leroy et al. (2008) < 
NGC 4449 is from Bolatto et al. (2008). Upper limits are at 5(7 signifi( 

in the mean because we are usually interested in how well 
a given number describes our whole sample, not how pre- 
cisely we have measured the mean. 

4. RESULTS 

Here we present our main observational results, how 
the star formation efficiency varies as a function of other 
quantities. We begin in §4.1 by showing the SFE as a 
function of three basic parameters: galactoccntric radius, 
stellar surface density, and gas surface density. Then in 
§4.2, we look at SFE as a function of the laws described 
in §2.1. Finally, in §4.3 we show the SFE as a function 
of the thresholds described in §2.2. 

We present these results as a scries of plots that each 
show SFE as a function of another quantity. These all 
follow the format seen in Figure 1, where we show SFE 
(y-axis) versus galactoccntric radius (x-axis), normalized 
to the optical radius, r25 . We plot the subsamples of spi- 
ral (top row) and dwarf galaxies (bottom row) separately. 

On the left, we show results for radial profiles. Each 
point shows the average SFE over one 10"-wide tilted 
ring in one galaxy. The color indicates whether the ISM 
averaged over the ring is mostly (> 50%) H I (blue) or 
H2 (magenta). Thick black crosses show all data binned 
into a single trend. For each bin, we plot the median, 
50% range (j/-error bar), and bin width (a;-error bar, here 
0.1 ras). 

On the right, we again show SFE as a function of ra- 
dius, this time calculated for each line of sight. We coadd 
all galaxies, giving equal weight to each, and pick con- 
tours that contain 90% (green), 75% (yellow). 50% (red), 
and 25% (purple) of the resulting data. Most numerical 
results use the annuli, which are easier to work with; 
these pixel-by-pixel plots verify that conclusions based 
on rings hold pixel-by-pixel down to kiloparsec scales. 

We do not analyze data with Egas < 1 M© pc~^ be- 



■ BIMA SONG (Heifer et al. 2003). NGC 3077 is from Walter et al. (2001), 
nice. 

cause the SFE is not well-determined for low gas surface 
densities; that is, we only address the question "where 
there is gas, is it good at forming stars?" Data with 
SsFR < 10~^ Mq yr^^ kpc^^ are treated as upper limits. 
These are red arrows in the radial profiles plots. In the 
pixel-by-pixel plots, hatched regions show the area inhab- 
ited by 95% of data with Esfr < 10"^ M© yr"! kpc^^ 
i.e., the hatched regions indicate the area where we are 
incomplete. In the pixel-by-pixel plots, we include data 
out to r25, while we plot radial profile data out to 1.2 r25. 



4.1. SFE and Other Basic Quantities 
4.1.1. SFE and Radius 

We argued in §1 that a critical observation for theories 
of galactic-scale star formation is that the SFE declines 
in the outer parts of spiral galaxies. Figure 1 shows this 
via plots of SFE against galactoccntric radius (normal- 
ized to ras) in our two subsamples. 

In spiral galaxies (top row), the SFE is nearly constant 
where the ISM is mostly H2 (magenta), which agrees with 
our observation of a linear relationship between Sh2 and 
EsFR in Bigiel et al. (2008). Typically, the ISM is equal 
parts H I and H2 at rgai = 0.43 ±0.18 r25 (§5.2). Outside 
this transition, the SFE decreases steadily with increas- 
ing radius. This decline continues to rgai > r25, the 
limit of our data. This is similar, though not identical, 
to the observation by Kennicutt (1989) and Martin & 
Kennicutt (2001) that star formation is not widespread 
beyond a certain radius. 

The SFE in spirals can be reasonably described in two 
ways. First, a constant SFE in the inner parts of galaxies 
followed by a break at 0.4 r25 (slightly inside the transi- 
tion to a mostly-H I ISM): 
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0.0 O.g 0.4 0.6 O.B 1.0 1.2 0.0 0.2 0.4 0.6 O.B 1.0 1.2 

Galactocentric Radius [r^^] Galactocenlric Radius [r^^] 



Fig. 1. — Star formation efficiency as a function of galactocentric radius in spiral (top row) and dwarf (bottom row) galaxies. The left 
panels show results for radial profiles; each point shows the average SFE over a 10"-wide tilted ring; magenta points are H2-dominated 
(Sh2 > Shi)i blue points are H I-dominated (Eh2 < Shi), and red arrows indicate upper limits. The right panels show data for individual 
lines of sight. We give each galaxy equal weight and choose contours that include 90%, 75%, 50%, and 25% of the data. The hatched 
regions indicate where we are incomplete. The top panels show a nearly fixed SFE in H2-dominated galaxy centers (magenta). Where 
H I dominates the ISM (blue), we observe the SFE to decline exponentially with radius; the thick dashed lines show fits of SFE to rgai 
(Equations 22 and 23). The vertical dotted line in the upper panels shows rgai at the H I-to-H2 transition in spirals, 0.43 ±0.18 r25 (§5.2). 



4.3 X 10-10 rgai < 0.4r25 _^ 

- j 2.2 X 10-^ exp (^) r,^, > OAr,, 

(21) 

Alternatively, we can adopt Equation 8, appropriate for 
a fixed SFE (H2), and derive the best-fit exponential re- 
lating i?niol to rgai, 

SFE - 5.25 X 10-10 TT^^ yr-i (22) 

^mol + 1 

-Rmoi = 10.6 exp (-rgai/0.21 r25) , 

which appears as a thick dashed line in the upper pan- 
els of Figure 1. The two ffis reproduce the observed 



SFE with similar accuracy; the scatter about each is 
« 0.26 dex, slightly better than a factor of 2. 

In dwarf galaxies (lower panels), we observe a steady 
decline in the SFE with increasing radius for all rgai, 
approximately described by 



SFE = 1.45 X 10-^ exp (-rgai/0.25 rzs) yr-^ (23) 

with ~ 0.4 dex scatter about the fit, i.e., a factor of 2-3. 

In dwarfs, we take Sgas ~ Shi, so that SFE= 
Esfr/Shi. For comparison with Equation 22, how- 
ever, we rewrite Equation 23 assuming that SFE (112)= 
5.25 X IQ-io yr-i, the value measured in spirals. In terms 
of i?moi, Equation 23 becomes 




Leroy 



optical Radius, [l<I>c] 

Fig. 2. — Stellar scale length, , as a function of isophotal radius, 
r25- Solid and dashed lines show r25 = (4.6 ± 0.8) U- 



SFE = = 5.25 X 10-10 R^oi yr"' (24) 

-Rinoi = 2.76 exp(-rgai/0.25 r25) • 

The outer parts of dwarfs, rgai > 0.4 r25, appear simi- 
lar to the outer disks of spiral galaxies in Figure 1. Sur- 
prisingly, however, we find the SFE to be higher in the 
central parts of dwarf galaxies than in the molecular gas 
of spirals. A higher SFE in dwarf galaxies is quite unex- 
pected. Their lower mctallicitics, more intense radiation 
fields, and weaker potential wells should make gas less 
efficient at forming stars. A simple explanation for the 
high observed SFE is the presence of a significant amount 
of H2. Figure 1 assumes that Sgas ~ Shi in dwarfs. If 
we miss a significant amount of H2 along a line of sight, 
we will overestimate the SFE because we underestimate 
Sgas- We quantify the possibility of substantial H2 in 
dwarfs in §5.3, but the magnitude of the effect can be 
read directly from Equation 24. At rgai = 0, if dwarf 
galaxies have the same SFE (H2) as spirals, Rmoi ~ 2.76, 
i.e., SH2 « 2.76 Shi. 

4.1.2. SFE and Stellar Surface Density 

Galactocentric radius is probably not intrinsically im- 
portant to a local process like star formation, but Fig- 
ure 1 suggests that local conditions covariant with radius 
have a large effect on the ability of gas to form stars. The 
radius, r25 , that we use to normalize the x-axis is defined 
by an optical isophote and thus measures stellar light. 
Therefore r25 is closely linked to the stellar distribution. 

Figure 2 shows this link directly. We plot stellar scale 
length, Z*, measured via an exponential fit to the 3.6/im 
profile as a function of r25 for our spiral subsamplc. We 
see that r25 — (4.6 ± 0.8) /, and that we could have 
cquivalently normalized the .x-axis in Figure 1 by Z*. We 
may suspect, then, that the stellar surface density, S*, 
underlies the well-defined relation between SFE and rgai 
observed in Figure 1. 

In Figure 3, we explore this connection by plotting SFE 
as a function of S*. In both spiral and dwarf galaxies, 
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we see a nearly linear relationship between SFE and S» 
where the ISM is H I-dominatcd (blue points). 

A basic result of THINGS is that over the optical disk 
of most star forming galaxies, the H I surface density 
varies remarkably little (Appendices E and F and Walter 
et al. 2008). Inspecting our atlas, one sees that Shi ~ 
6 Mq pc^ (within a factor of 2) over a huge range of 
local conditions, including most of the optical disk in 
most galaxies. Because Sgas is nearly constant in the 
H I-dominated (blue) regime, SFE oc S* approximately 
defines a line of fixed specific star formation rate (SSFR), 
i.e., star formation rate per unit stellar mass. 

The inverse of the SSFR is the stellar assembly time, 
T* = S^/SgiTfl. This is the time required for the present 
star formation rate to build up the observed stellar disk. 
In our spiral subsample, the mean log^o t* « 10.5 ± 0.3, 
i.e., 3.2 X 10^" years or slightly more than 2 Hubble times. 
Dwarf galaxies have shorter assembly times, log^Q r* w 
10.2 ± 0.3 years, about a Hubble time (dashed lines in 
Figure 3 show these values using average values of Sgas 
for each subsample). Taking these numbers at face value, 
dwarfs are forming stars at about their time-average rate, 
while spirals arc presently forming stars at just under half 
of their average rate. 

We only observe SFE oc S* where the ISM is mostly 
H I. Where the ISM is mostly H2 in spirals galaxies, 
we observe a constant SFE at a range of S,; similar to 
the constancy as a function of rgai observed in the inner 
parts of spirals (Figure 1). The transition between these 
two regimes occurs at S* = 81 ± 25 Mq pc~^ (§5.2) 
in spirals. In dwarfs, lines of sight with S* above this 
transition value exhibit systematically high SFE, lending 
further, albeit indirect, support to the idea that these 
points correspond to unmeasured H2. 

Figures 1 and 3 show that where the ISM is mostly 
H2, the star formation rate per unit gas (SFE) is nearly 
constant and that where the ISM is mostly H I, the star 
formation rate per unit stellar mass (SSFR) is nearly 
constant. Together these observations suggest that H2, 
stars, and star formation have similar structure with all 
three embedded in a relatively flat distribution of H I. 
Figure 4 shows that the scale lengths of these three dis- 
tributions are, in fact, comparable. The star formation 
rate (black) and CO (gray) scale lengths of spiral galaxies 
are both roughly equal to the stellar scale length: 

Ico = (0.9 ± 0.2) h and Zsfr = (1 ± 0.2) h ■ (25) 

Regan et al. (2001) also found that Zco ~ I* comparing 
if -band maps to BIMA SONG and Young et al. (1995) 
found Ico ~ 0-2 r25, which is almost identical to our 
Ico ~ 0.9 h and (4.6 ± 0.8) h = r25- 

4.1.3. SFE and Gas Surface Density 

This link between S* and the SFE is somewhat sur- 
prising because it is common to view Ssfr, and thus the 
SFE, as set largely by Sgas alone over much of the disk of 
a galaxy (following, e.g., Kennicutt 1998a). In Figure 5 
we show this last slice through SFR-stars-gas parameter 
space, plotting SFE as a function of Sgas- 

As in Figures 1 and 3, wc observe two distinct regimes. 
In spirals, where Sgas > 14 ± 6 Mq pc^^ (§5.2) the 
ISM is mostly H2 and we observe a fixed SFE. This 
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Fig. 3. — SFE as a function of stellar surface density, S,, in spiral (top row) and dwarf (bottom row) galaxies. Conventions and symbols 
are as in Figure 1. Dashed diagonal lines show the linear relationship between SFE and S, expected for the mean stellar assembly time and 
Sgaa for each subsample. Vertical dotted lines show S* where the ISM is equal parts H I and H2 in spirals (§5.2), S« = 81 ± 25 Mq pc~^. 



Sgas, shown by a vertical dotted line, corresponds ap- 
proximately to both N{H) ^ 10^^ cm~^ star forma- 
tion threshold noted by Skillman (1987) and the satura- 
tion value for H I observed by, e.g., Martin & Kennicutt 
(2001) and Wong & Blitz (2002) (and seen strikingly in 
THINGS at Sgas = 12 Mg pc"^ by Bigiel et al. 2008, who 
quote Egas = 9 M© pc~^ but do not include helium). 

In contrast to rgai and S*, Egas does not exhibit a clear 
correlation with the SFE where the ISM is mostly H I. 
Instead, over the narrow range Egas ~ 5-10 Mq pc~^, 
the SFE varies from ~ 3 x 10"" to 10"^ yr"!. We see 
little evidence that Ehi plays a central role regulating the 
SFE in either spirals or dwarfs. Rather, the most striking 
observation in Figure 5 is that Ehi exhibits a narrow 
range of values over the optical disk and is therefore itself 
likely subject to some kind of regulation. 

The possibility of a missed reservoir of molecular gas 
in dwarfs is again evident from the lower panels in Figure 



5. A subset of data has SFE higher than that observed 
for II2 in spirals and just to the left of the H I saturation 
value. If H2 were added to these points, they would move 
down (as the SFE decreases) and to the right (as Egas 
increases), potentially yielding a data distribution similar 
to that we observe in spirals. 

4.2. SFE and Star Formation Laws 

We now ask whether the star formation laws proposed 
in §2.1 can explain the radial decline in SFE and whether 
SFE (II2), already observed to be constant as a function 
of Tgai, E,, and Egas (but with some scatter), exhibits 
any kind of systematic behavior. We compare the SFE to 
four quantities that drive the predictions in Table 1: gas 
surface density (already seen in Figure 5), gas pressure 
(density) , the orbital timescale, and the derivative of the 
rotation curve, f3. 

4.2.1. Free-Fall Time in a Fixed Scale Height Disk 
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Fig. 4. — The scale lengths of star formation (blaek) and CO 
(gray) as a funetion of the stellar scale length (x-axis). All three 
scale lengths are similar, the dashed lines show slope unity and 
±30% (the approximate scatter in the data). 

A dashed line in Figure 5 illustrates SFE oc Sgas, ex- 
pected if the SFE is proportional to the free-fall time in 
a fixed scale height gas disk (similar to the Kcnnicutt- 
Schmidt law, Kennicutt 1998a). The normalization 
matches the H2-dominated parts of spirals and roughly 
bisects the range of SFE observed for dwarfs, but large 
areas of the disk have much lower SFE than one would 
predict from this relation. Adjusting the normalization 
can move the line up or down but cannot reproduce the 
distribution of data observed in Figure 5. 

The culprit here is the small dynamic range in Ehi- 
Because Shi does not vary much across the disk, while 
the SFE docs, the free fall time in a fixed scale height 
disk, or any other weak dependence of SFE on Sgas alone, 
cannot reproduce variations in the SFE where the ISM 
is mostly H I. A quantity other than Egas must play an 
important role at radii as low as 0.5 r25 (a fact already 
recognized by Kennicutt 1989, among others). 

4.2.2. Free-Fall Time in a Variable Scale Height Gas 
Disk; Pressure and ISM Phase 

We saw in §4.1 that where the ISM is mostly H I, the 
SFE correlates better with than with Sgas- This 
might be expected if the stellar potential well plays a 
central role in setting the volume density of the gas, pgas, 
because varies mtich more strongly with radius than 
Ehi- In §2 we present two predictions relating SFE to 
Pgas: that the timescale over which GMCs form depends 
on the Tfi, the free-fall time in a gas disk with a scale 
height set by hydrostatic equilibrium^, and that the ra- 
tio i?nioi = 5]h2/^hi depends primarily on midplane gas 
pressure, Ph- 

Under our assumption of a fixed (Tgas, Ph oc pgas and 
both predictions can be written as a power law relating 
SFE or iimoi to Ph- In Figure 6 we plot SFE as a function 

of pgas and Ph (top and bottom .T-axis) , estimated from 
hydrostatic equilibrium (Equation 9). 

® Hereafter Tff refers only to the free fall time in a gas disk with 
a scale height set by hydrostatic equilibrium. 



Where the ISM is mostly H2 (magenta points) in spi- 
rals (top row), we observe no clear relationship between 
Ph and SFE, further evidence that SFE (H2) is largely 
decoupled from global conditions of the ISM in our data. 

Where the ISM is mostly H I (blue points) in dwarf 
galaxies and the outer parts of spirals, the SFE correlates 
with Ph. Ph predicts the SFE notably better than Egas in 
this regime, supporting the idea that the volume density 
of gas (at least H I) is more relevant to star formation 
than surface density. Wong & Blitz (2002) and Blitz 
& Rosolowsky (2006) observed a continuous relationship 
between Pmoi and Ph, mostly where Eh2 ^ Ehi- Figure 
6 suggests that such a relationship extends well into the 
regime where H I dominates the ISM. 

The solid line in Figure 6 illustrates the case of 1% of 
the gas formed into stars per Tff (SFE cx Pgas), a typ- 
ical value at the H I to II2 transition in spirals (§5.2). 
Adjusting the normalization slightly, such a line can in- 
tersect both the high and low end of the observed SFE in 
spirals, but predicts variations in SFE (H2) that we do 
not observe and is too shallow to describe dwarf galaxies. 

The dash dotted line shows P,noi fx t^^ oc P]'^"'', ex- 
pected for GMC formation over a free fall time. In dwarf 
galaxies, where we take Egas = Shi, this is equivalent 
too SFEoc ■ This description can describe spirals at 
high and intermediate Ph, but is too shallow to capture 
the drop in SFE at large radii in spirals and across dwarf 
galaxies. If ts is the characteristic timescale for GMC 
formation, effects other than just an increasing timescale 
must suppress cloud formation in these regimes. 

A dashed line shows the steeper dependence, Pmoi oc 
PjJ'^, expected for low Pmoi based on modeling by 
Elmegreen (1993). This may be a reasonable description 
of both spiral and dwarf galaxies (note that at high SFE, 
Ph may be underestimated in dwarf galaxies because we 
fail to account for H2). We explore how Ph relates to 
Rmoi more in §5. 

4.2.3. Orbital Timescale 

The orbital timescale, Torb, varies strongly with radius 
and Kennicutt (1998a) found Torb to be a good predictor 
of disk-averaged SFE. In Figure 7, we plot SFE as a 
function of Torb in our sample. 

The solid line shows 6% of the gas converted to stars 
per Torb and is a reasonable match to spirals near the 
H I-to-H2 transition (vertical dotted line) . This value 
agrees with the range of efficiencies found by Wong & 
Blitz (2002) and with Kennicutt (1998a), who found w 
7% of gas converted to stars per Torb averaged over galaxy 
disks (converted to our adopted IMF). Like Wong & Blitz 
(2002), we do not observe a clear correlation between 
SFE and Torb where the ISM is mostly H2. 

Where the ISM is mostly HI (blue points), the SFE 
clearly anti-correlates with Torb in both spiral and dwarf 
galaxies. However, we do not observe a constant effi- 
ciency per Torb- In both subsamplcs, SFE drops faster 
than Torb increases, so that data at large radii (longer 
''orb, lower SFE) show lower efficiency per Torb than those 
from inner galaxies. Although Torb correlates with the 
SFE, the drop in Torb is not enough on its own to explain 
the drop in SFE. 

We reach the same conclusion if we posit that Torb 
is the relevant timescale for GMC formation, so that 
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Fig. 5. — SFE as a function of Egas in spiral (top row) and dwarf (bottom row) galaxies. Conventions and symbols are the same as in 



Figure 1. The vertical dotted line shows Sgas at the H I-to-H2 transition in spirals (§5.2), Eg 



14±6 M, 



pc 



^ . The dashed line shows 



the SFE proportional to the free-fall time in a fixed scale height disk. Clearly the line cannot describe both high and low SFE data, even 
if the normalization is adjusted, and so changes in this timescale cannot drive the radial decline that we observe in the SFE. 



^moi oc T~j^. The dashed hnes in Figure 7 show this 
relation combined with a fixed SFE (H2) and normalized 
to i?i„oi = 1 at Torb = (1-8 ± 0.4) X 10^ years, which 
we observe at the H I-to-H2 transition in spirals (§5.2). 
This dependence is even shallower than SFE oc t^^^ and 
cannot reproduce the SFE in both inner and outer disks 
by itself. If Torb is the relevant timescale for cloud for- 
mation, then the fraction of gas that is actively forming 
stars must vary substantially between the middle and the 
edge of the optical disk. 

4.2.4. Derivative of the Rotation Curve, (3 

Tan (2000) suggests that cloud-cloud collisions regu- 
late the SFE. The characteristic timescale for such colli- 
sions is Torb modified by the effects of galactic shear. We 
saw in Figure 7 that the SFE of molecular gas is not a 
strong function of Torb- Therefore, in Figure 8, we plot 
the SFE as a function of /3, the logarithmic derivative of 



the rotation curve (we plot SFE against Qgas, the other 
component of this timescale in §4.3.1). This isolates the 
effect of differential rotation; /3 = for a flat rotation 
curve and (3 = 1 for solid body rotation (no shear). 

Figure 8 shows a simple relationship between /3 and 
SFE in spirals: /3 > is associated with high SFE. High 
(3 occurs almost exclusively at low radius (where the ro- 
tation curve rises steeply) and in these regions the ISM 
is mostly H2 with accordingly high SFE. On the other 
hand, the outer disks of spirals have /3 ^ and a wide 
range of SFE. Beyond basic relationship, it is unclear 
that (3 has utility predicting the SFE. In particular, we 
see no clear relationship between SFE and (3 where the 
ISM is mostly H2 (magenta points). If collisions between 
bound clouds regulate the SFE, we would expect an anti- 
correlation between (3 and SFE because cloud collisions 
are more frequent in the presence of greater shear. 

In dwarf galaxies increasing (3 corresponds mostly to 
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Fig. 6. — SFE as a function of midplane hydrostatic gas pressure, Ph (bottom x-axis) and equivalent volume density (top a;-axis) in 
spiral (top row) and dwarf (bottom row) galaxies. Conventions follow Figure 1. The vertical line shows Ph at the H I-to-H2 transition 
in spirals, logjg Ph/^B [K cm~^] 4.36. The solid line illustrates 1% of gas converted to stars per disk free-fall time. Dash-dotted and 
dashed lines show -Rmol oc P^'^ (''^ff ^) ^^'^ ^mol °^ (Elmegreen 1993). For our adopted o-gas = 11 km s~^ and including helium: 

p [g cm-3] = 1.14 X 10~2® (Ph/kB) [K cm-3] and n [cm'^] = 4,4 x 10^3 p [g cm-^] . 



increasing SFE. This relationship has the sense of the 
shear threshold proposed by Hunter et al. (1998a), that 
where rotation curves are nearly solid body low shear al- 
lows clouds to form via instabilities aided by magnetic 
fields (see also Kim & Ostriker 2001). The rotation 
curves in dwarf galaxies rise more slowly than those in 
spirals, leading to /3 > over a larger range of radii in 
dwarf galaxies and limiting = to the relative out- 
skirts of the galaxy. A positive correlation between /? 
and SFE is opposite the sense expected if cloud colli- 
sions are important: at high /3 collisions should be less 
frequent. 



4.3. SFE and Thresholds 



The decline in the SFE where the ISM is mostly H I is 
too dramatic to be reproduced across our whole sample 
by changes in Torb or Tg alone. This may be because at 
large radii a significant amount of gas is simply unrelated 
to star formation. If the fraction of gas that is unable to 
form GMCs increases with radius, the SFE will decline 
independent of any change in GMC formation time. Here 
we consider the SFE as a function of proposed star for- 
mation thresholds: gravitational instability in the gas 
alone (Qgas), in a disk of gas and stars (Qstars+gas), the 
ability of instabilities to develop before shear destroys 
them, and the ability of a cold gas phase to form. 

First we plot each threshold as a function of galac- 
tocentric radius in spiral (Figure 9) and dwarf galax- 
ies (Figure 10). Individual points correspond to aver- 
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Fig. 7. — SFE as a function of the orbital timescale, Torbi in spiral (top row) and dwarf (bottom row) galaxies, following the conventions 
from Figure 1. The solid line shows 6% of gas converted into stars per Torb- The dashed line shows the expected SFE if -Rmol = Sh2/Shi oc 
T~}. The SFE is a well-defined function of To^bi but the decline in t^j-i, alone cannot reproduce the radial decline in SFE or -Rmol- 



ages over 10"-wide tilted rings. For magenta points 
Eh2 > Shi and for blue points Sh2 < Shi- The gray 
region in each plot shows the nominal condition for in- 
stability, i.e., where we expect star formation to occur. 
Red arrows indicate data outside the range of the plot. 

We proceed creating plots like Figure 1 for each thresh- 
old and comparing them to Figure 9. We expect super- 
critical gas to exhibit a (dramatically) higher SFE than 
subcritical gas, where star formation proceeds only in 
isolated pockets or not at all. 

4.3.1. Gravitational Instability in the Gas Disk 

Figure 11 shows SFE as a function of Qgas, Toomre's 
Q parameter for a thin gas disk; the top left panels in 
Figures 9 and 10 show Qgas as a function of radius. 

In each plot, a gray area indicates the theoretical con- 
dition for instability. We see immediately that almost 
no area in our sample is formally unstable. Rather, most 
lines of sight are strikingly stable, Qgas '--^ 4 is typical 



inside ^ 0.8 and Qgas > 10 is common. 

We find no clear evidence for a Qgas threshold (at any 
value) that can unambiguously distinguish regions with 
high SFE from those with low SFE. In spirals, Qg^s ^ 2.5 
appears to be a sufficient, but by no means necessary 
condition for high SFE; there are also areas where the 
ISM is mostly H2, SFE is quite high and Qgas > 10. 
In dwarfs Qgas appears, if anything, anti-correlated with 
SFE, though this may partially result from incomplete 
estimates of Sgas. 

These conclusions appear to contradict the findings by 
Kcnnicutt (1989) and Martin & Kennicutt (2001), who 
found marginally stable gas (Qgas 1-5) across the op- 
tical disk with a rise in Qgas corresponding to dropping 
SFE at large radii. In fact, after correcting for differ- 
ent assumptions, our median Qgas matches theirs quite 
well. Both Kennicutt (1989) and Martin & Kennicutt 
(2001) assumed Xco = 2.8 x 10^° cm^^ (K km s"^)"^ 
and (Tgas = 6 km s~^, while we take Xco — 2.0 x 
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Fig. 8. — SFE as a function of /3, the logarithmic derivative of the rotation curve in spiral (top row) and dwarf (bottom row) galaxies. 
/3 = 1 for solid body rotation and /3 = for a fiat rotation curve. If collisions between GMCs were important to triggering star formation, 
we would expect the SFE in the H2 dominated (magenta) parts of spirals to be higher for low f3 (high shear), which is not apparent from 
the data. 



10^° cm~^ (K km s~-^)~^ and Cgas = H km s~^. As a re- 
sult, we estimate less H2 and more kinetic support than 
they do for the same observations. If we match their as- 
sumptions, our median Qgas in spirals and the outer parts 
of dwarfs agrees quite well with their threshold value, 
though we find the central regions of dwarfs systemati- 
cally above this value (as did Hunter et al. 1998a). We 
show this in Figures 9, 10, and 11 by plotting the Mar- 
tin & Kennicutt threshold converted to our assumptions 
(Qgas ^ 3.9) as a a dashed line. 

The main observational difference between our result 
and Martin & Kennicutt (2001) is that Qgas shows much 
more scatter in our analysis. As a result, a systematic 
transition from low to high Qgas near the edge of the 
optical disk is not a universal feature of our data, though 
a subset of spiral galaxies do show increasing Qgas at 
large radii (Figure 9). 

This discrepancy in Qgas derived from similar data 



highlights the importance of assumptions. The largest 
effect comes from CTgas, which we measure to be ~ 
11 km s~^ and roughly constant in H I~dominated outer 
disks (Appendix B). We assume Cgas to be constant ev- 
erywhere, an assumption that may break down on small 
scales and in the molecular ISM. In this case we expect 
(Tgas to be locally lower than the average value, lower- 
ing Qgas and making gas less stable. Black dots in the 
upper right panel of Figure 11 show the effect of chang- 
ing CTgas from 11 km s~^ (our value) to 6 km s^^ (the 
Martin & Kennicutt value) and then to 3 km s~^, the 
value expected and observed for a cold H I component 
(e.g. Young et al. 2003; Schaye 2004; de Blok & Walter 
2006). If most gas is cold then Qgas may easily be < 1 
for this component (if only a small fraction of gas is cold, 
the situation is less clear). 

4.3.2. Gravitational Instability Including Stars 
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Fig. 9. — Radial behavior of thresholds in spiral galaxies: {top left) gravitational instability due to gas self gravity; {top right) gravitational 
instability due to the combination of self-gravity and stellar gravity; {bottom left) competition between cloud formation and destruction by 
shear; {bottom right) formation of a cold phase. Each point shows average Scrit/Sgas over one 10" tilted ring in one galaxy. In magenta 
rings, the ISM is mostly H2, in blue rings the ISM is mostly H I. Gray regions show the condition required for star formation. 



Stars dominate the baryon mass budget over most of 
the areas we study and stellar gravity may be expected 
to affect the stability of the gas disk. In §2 we described 
a straightforward extension of Qgas to the case of a disk 
containing gas and stars (Rafikov 2001). In Figure 12 
we plot SFE as a function of this parameter, Qstars+gas , 
which we plot as a function of radius in the top right 
panels of Figures 9 and 10. 

The gray region indicates where gas is unstable to ax- 
isymmetric collapse. Including stars does not render 
large areas of our sample unstable, but it does imply that 
most regions are only marginally stable, Qstars+gas ~ 1-6. 
This in turn suggests that it is not so daunting to induce 
collapse as one would infer from only Qgas- 

In addition to lower values, Qstars+gas exhibits a much 
narrower range of values than Qgas, mostly areas in both 
spiral and dwarf galaxies show Qstars+gas — 1-3 - 2.5. 
This may offer support to the idea of self-regulated star 



formation, but it also means that Qstars+gas offers little 
leverage to predict the SFE. High SFE, mostly molecular 
regions show the same Qstars+gas as low SFE regions from 
outer disks (indeed, the highest values we observe come 
from the central parts of spiral galaxies). 

As with Qgas, our assumptions have a large impact on 
Qstars+gas- In addition to Cgas and Xco (which affect the 
calculation via Qgas), the stellar velocity dispersion, cr*, 
and mass-to-light ratio, T^, strongly affect our stability 
estimate. We assume that cr* oc E°'^ in order to yield a 
constant stellar scale height. If we instead fixed a^,, we 
would derive Qstars+gas increasing steadily with radius. 
Radial variations in may create a similar effect. 

Boissier et al. (2003) find similar results to our own 
when they incorporate stars in their stability analysis; 
they adopt a lower cTgas than we do, but also lower Xqo 
and the effects roughly offset. Yang et al. (2007) recently 
derived Qstars+gas across the LMC and found widespread 
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Fig. 10. — Radial behavior of thresholds in dwarf galaxies: (top left) gravitational instability due to gas self gravity; (top right) 
gravitational instability due to the combination of self-gravity and stellar gravity; (bottom left) competition between cloud formation and 
destruction by shear; (bottom right) formation of a cold phase. Each point shows average Scrit/Sgas over one 10" tilted ring in one galaxy. 
Gray regions show the condition required for star formation. 



instability that corresponded well with the distribution of 
star formation. If we match their adopted agas (5 km s^^) 
and assumptions regarding cr, (constant at 15 km s~^) 
we also find widespread instability throughout our dwarf 
subsample, Qstars+gas decreasing with radius; we find a 
similar result for spirals if we fix a typical outer-disk ct» . 

Our approach is motivated by observations of disk 
galaxies (see Appendix B), but direct observations of cr, 
at large radii are still sorely needed. 

4.3.3. Shear Threshold 

If clouds form efficiently, e.g., through the aid of mag- 
netic fields to dissipate angular momentum, then Hunter 
et al. (1998a) suggest that the time available for a pertur- 
bation to grow in the presence of destructive shear may 
limit where star formation is widespread. Kim & Os- 
triker (2001) describe a similar scenario where magneto- 
Jeans instabilities can grow in regions with weak shear or 



strong magnetic fields. In the bottom left panels of Fig- 
ures 9 and 10, we plot this shear threshold as a function 
of radius and in Figure 13 we compare it to the SFE. 

The gray region shows the condition for instabilities 
to grow into GMCs, Scrit,A/Sgas < 1- This matches the 
condition Qgas < 1 where /3 = 0, e.g., in outer disks of 
spirals. In the inner parts of spirals and in dwarf galax- 
ies, however, Ecrit.A/Sgas is lower than Qgas, i-e., the 
conditions for star formation are more nearly supercriti- 
cal (because shear is low in these regions). These areas 
harbor II2 or widespread star formation, so supercritical 
values are expected. 

This trend of more supercritical data at lower radii 
agrees with the steady increase of SFE with decreasing 
radius that we saw in Figure 1. However, the scatter 
in Scrit.A/Sgas < 1 is as large as that in Qgas (as one 
would expect from their forms, see Table 1). As a result, 
a direct plot of SFE against Ecrit,A/Sgas does not yield 
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Fig. 11. — SFE as a function of Qgas, the Toomre Q parameter, which measures instability to axisymmetric collapse in a gas disk. 
Symbols and conventions follow Figure 1. The gray region shows where instability is expected. A dashed line in the top right panel shows 
the Qgas threshold derived from Ha emission by Martin & Kennicutt (2001) converted to our assumptions. In the same panel, we show 
the effect on Qgas of changing erg from our adopted 11 km s~^ to 6 km s~^ and then 3 km s~^, expected for a cold phase. 



clear threshold behavior or a strong correlation between 
Scrit,A/Sgas and SFE. The strongest conclusion we can 
draw is that the inner parts of both spiral and dwarf 
galaxies are marginally stable for the shear threshold (an 
improvement over Qgas and Qstars+gas in these regions). 

4.3.4. Cold Phase Formation 

Even where the ISM is stable against gravitational col- 
lapse on large scales, star formation may still proceed if 
a cold (narrow-line width) phase can form locally and 
thus induce gravitational instability in a fraction of the 
gas (recall the effect of lower CTgas in Figure 11). Schaye 
(2004) argued that this is the usual path to star for- 
mation in the outer parts of galaxies and modeled the 
critical gas surface density for such a phase to form, 
Eso4- The bottom right panels in Figures 9 and 10 show 
Eso4/5jgas as a function of radius and Figure 14 shows 
the SFE as a function of this ratio. The gray area in 



both figures shows where a cold phase can form. 

We calculate Sso4 from Equation 20, which depends on 
I/[10^ cm~^ s~^], the flux of ionizing photons. In outer 
disks, we assume / = 10^ cm^^ s^^, Schaye's fiducial 
value, and in inner disks we take / cx Esfr, 

/ « 10^ cm-2 s-i ( ^) . (26) 

V5 X 104 M© yr-i kpc"V 
The normalization is the average Ssfr between 0.8- 
1.0 in our spiral subsample. 

Equation 20 also accounts for variations about Schaye's 
fiducial metallicity Z = 0.1 Zq, typical for the outer disk 
of a spiral. We lack estimates of Z and so neglect this 
term but note the sense of the uncertainty. Inner galaxy 
disks will tend to have higher metallicities, which will 
lower Sso4- We already find Sgas > '^soi over most inner 
disks; therefore missing Z seems unlikely to seriously bias 
our results. 
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Fig. 12. — SFE as a function of Qstars+gas (Rafikov 2001), which measures instability in a gas disk in the presence of a coUisionless 
stellar disk. Symbols and conventions follow Figure 1. The gray region indicates where gas is unstable. Compared to Qgas, including stars 
renders the disk more nearly unstable and yields a much lower range of values. 



Figures 9,10, and 14 show that we expect a cold phase 
over most of the disk in both spiral and dwarf galaxies. 
In our spiral subsample, most data inside rgai ~ 0.9 
meet this criterion. Because most subcritical data come 
from large radii, we also find that most lines of sight with 
Sgas < SS04 exhibit low SFEs or upper limits. 

Because most data are supercritical, the Schaye (2004) 
threshold is of limited utility for predicting the SFE 
within a galaxy disk. Schaye (2004) does not predict the 
ratio of H2-to-H I where cold gas forms; he is primar- 
ily concerned with the edges of galaxies. Figures 9 and 
14 broadly confirm that his proposed threshold matches 
both the edge of the optical disk and the typical thresh- 
old found by Martin & Kennicutt (2001). 

This relevance of this comparison to the SFE within 
the optical disk is that based on the Schaye (2004) model, 
we expect a widespread narrow-line phase throughout 
most of our galaxies (Wolfire et al. 2003, obtain a similar 
result for the Milky Way) . This suggests that cold phase 



formation followed by collapse may be a common path 
to star formation and offers a way to form stars in our 
otherwise stable disks. 

5. DISCUSSION 

In §4.1, 4.2 and 4.3 we examined the SFE as a function 
of basic physical parameters, laws, and thresholds. Here 
we collect these results into general conclusions regard- 
ing the SFE in galaxies and identify key elements of a 
successful theory of star formation in galaxies. 

5.1. Fixed SFE 0/H2 

Using a data set that overlaps the one presented here, 
Bigiel et al. (2008) found a linear relationship between 
EsFR and SH2- Here we extend that finding: Where the 
ISM is mostly H2 in spiral galaxies, the SFE does not 
vary strongly with any of the quantities that we consider, 
including radius, Sgas, S,, Ph, T-orb, and /3. We plot 
SFE (H2) as a function of each of these quantities in 
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Fig. 13. — SFE as a function of Scrit,A/Sgas, the threshold for cloud growth in the presence of shear (Hunter et al. 1998a) for spiral 
(top row) and dwarf (bottom row) galaxies. Conventions and symbols follow Figure 1. The gray area shows where clouds should be able 
to survive distribution by shear. 



Figure 15. The median value for tilted rings from our 
spiral subsample is log-^„ SFE (H2) = -9.28 ± 0.17, i.e., 

SFE (H2) = 5.25 ±2.5 X 10^1° yr~i . (27) 

Constant SFE (H2) might be expected if 1) conditions 
within a GMC, rather than the larger scale properties 
of the ISM, drive star formation (e.g., Krumholz & Mc- 
Kee 2005) and 2) GMC properties are relatively univer- 
sal rather than, e.g., a sensitive function of formation 
mechanism or environment. This appears to be the case 
in the inner Milky Way (excluding the Galactic center) 
and in M31 and M33, where GMC properties are largely 
a function of cloud mass alone (Solomon et al. 1987; 
Rosolowsky et al. 2003; Rosolowsky 2007; Blitz et al. 
2007; Bolatto et al. 2008). The constancy of SFE (H2) 
hints that a similar case holds in our spiral subsample. 

Figure 6 illustrates why (relatively) universal GMC 
properties may be plausible in our sample. From Equa- 
tion 9, the internal pressure of a starless GMC with 



Sgas ~ 170 Mq pc-2 (Solomon et al. 1987) is Ph/fe 
10^ K cm""^. This is the highest value we plot in Figures 
6 and 15 and only a small fraction of our data have higher 
Ph so that even where the ISM is mostly H2, Ph is usu- 
ally well below the typical internal pressure of a GMC. 
Thus GMCs are not necessarily pressure-confined, which 
allows the possibility of bound, isolated GMCs out of 
pressure equilibrium with the rest of the ISM. In this 
case, the environmental factors that we consider may 
never be communicated to GMCs (though some mech- 
anism may still be needed to damp out any imprint left 
by environment during GMC formation). 

The range of Ph in our sample also underscores that 
one should not expect a constant SFE (H2) to extend 
to starburst conditions, where Ph and Sgas on kiloparsec 
scales exceed those found for individual Galactic GMCs 
and SFE (H2) is observed to vary strongly with local 
conditions (e.g., Kennicutt 1998a; Riechers et al. 2007). 

Another important caveat is that the distribution 
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Fig. 14. — SFE as a function of Eso4/Egas, the threshold for the formation of a cold phase Schaye (2004), for spiral (top row) and dwarf 
galaxies (bottom row). Conventions and symbols follow Figure 1. The gray area indicates where Schaye (2004) estimates that a cold phase 
can form. Most areas where we observe star formation meet this criterion and the areas that do not tend to have low SFE. 



of GMC masses is observed to vary with environment 
(Rosolowsky 2005), possibly as a result of varying for- 
mation mechanisms. This suggests that either the SFE 
of a GMC is only a weak function of its mass (and thus 
other properties) or that real variations in SFE (H2) may 
exist in dwarf galaxies and the outskirts of spirals. 

5.2. Conditions at the H I-to-H2 Transition in Spirals 

In spiral galaxies, the transition between an H I- 
dominated ISM and a mostly-Il2 ISM occurs at a char- 
acteristic value for most quantities. This can be seen 
from Figures 1, 3, 5, 6, and 7, in which H I-dominated 
regions (blue points) typically occupy one region and H2- 
dominated regions (magenta points) occupy another. 

Table 5 gives our estimates of properties where Shi ~ 
I]h2 in spiral galaxies. For each galaxy, we measure the 
median of the property in question over all pixels where 
Sh2 = 0.8 - 1.2 Ehi- Table 5 lists the median transition 
value in our spiral subsample, along with the (la) scatter 



and log scatter among galaxies. These values appear as 
dotted vertical lines in Figures 1, 5, 3, 6, and 7. Note 
that methodology — the choice to use pixels or rings, to 
interpolate, use the mean or median, etc. — affects the 
values in Table 5 by - 20%. 

From Table 5, we find that physical conditions at the 
H I-to-H2 transition are fairly similar to those found 
in the solar neighborhood. The orbital time is « 
1.8 X 10^ years and the free-fall time in the gas disk 
is w 4.2 X 10^ years. The midplane gas pressure is 
Ph/fce ~ 2.3 X IC* cm^'^ K, corresponding to a par- 
ticle density n ~ 1 cm^"^. The baryon mass budget 
in the disk is dominated by stars, « 81 Mq pc~^ 
while Egas ~ 14 M0 pc~^. Accordingly, the gas is 
stable against large scale gravitational collapse on its 
own (Qgas ~ 3.8), but in the presence of stars is only 
marginally stable Qstars+gas ~ 1-6. 

Approximately 1% of gas is converted to stars per free 
fall time at the transition, in agreement with expecta- 
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Fig. 15. — SFE (H2) in individual tilted rings from spiral galaxies as a function of: {top left) galactocentric radius, (top middle) H2 surface 
density, {top right) stellar surface density, {bottom left) midplane pressure, {bottom middle) orbital timescale, (bottom right) logarithmic 
derivative of the rotation curve. Gray lines show the median logj^Q SFE (H2) = —9.28 ± 0.17 for our data. 

in dwarfs indeed exceeds SFE (H2) in spirals. We argue 
that an unaccounted-for reservoir of H2 is a more likely 
explanation. 

The SFE (H2 ) that we observe in spiral galaxies offers 
an approximate way to estimate how much H2 may be 
present. If we assume that SFE (H2) is the same in dwarf 
and spiral galaxies then we can calculate Eh2 from the 
observed Ssfr via 



TABLE 5 

Conditions at the H I-T0-H2 Transition 



Quantity 




Median 


Scatter 


Scatter 






Value'' 




in logio 


r-gal [?-25] 




0.43 


0.18 


0.17 


S. [Mq pc-2 




81 


25 


0.15 


Sgas [Mq pc" 




14 


6 


0.18 


Ph/kB [cm-3 


K] 


2.3 X lO'' 


1.5 X 10"* 


0.26 


m [yr] 




4.2 X lO'^ 


1.2 X lO'^ 


0.14 


T-orb [yr] 




1.8 X 10* 


0.4 X 10* 


0.09 


Qgas 




3.8 


2.6 


0.31 


Qstars+gas 




1.6 


0.4 


0.09 



''Median value in the spiral subsample. 

tions by Krumholz & McKee (2005). About 6% of gas 
is converted to stars per Torb- This agrees well with 
the disk-averaged value of ~ 7% derived by Kennicutt 
(1998a) (adapted to our IMF and C0-to-H2 conversion 
factor) and with the range of efficiencies found by Wong 
& Blitz (2002). 

5.3. H2 in Dwarf Galaxies 

Because of uncertainties in Xco^ we do not directly 
estimate the amount of H2 in dwarf galaxies . However, 
indirect evidence suggests that a significant part of the 
ISM is H2 in the central parts of these galaxies. Specif- 
ically, we observe very high SFE in the centers of dwarf 
galaxies — higher than SFE (H2) in spirals — often un- 
der conditions associated with an H2-dominated ISM in 
spirals (§5.2). It would be surprising if the SFE of H I 



10' 



^H2 



■^SFR 



5.25 X 10 



-10 



(28) 



This treatment suggests that in our typical dwarf 
galaxies, most of the ISM is H2 within 0.25 r25- This 
may be seen directly from Equation 24, which translates 
our fit of SFE to radius to a relation between i?moi and ra- 
dius assuming Equation 28. From Equation 24, Sh2/5]hi 
in dwarf galaxies is 1-2 inside ~ 0.25 r25, rising as high 
as ~ 3 at Tgai = 0. 

5.4. Environment- Dependent GMC/ii2 Formation 

Where the ISM is H I-dominated — in dwarf galaxies 
and outside the H I-to-Il2 transition in spirals — the SFE 
declines steadily with increasing radius. In this regime, 
the SFE is covariant with a number of environmental 
factors, including S*, pressure, density, free fall time, 
and orbital timescale. This observation, together with 
those in §5.1 and 5.2, implies that while star formation 
within GMCs is largely decoupled from environment, the 
formation of H2 / GMCs from H I depends sensitively on 
local conditions. 
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Fig. 16. — The H2-to-H I ratio, /?niol = Sh2/Shi, as a function of {top left) radius, (top right) E,, (bottom left) Pjj (oc Tg ), and (bottom 
right) r^rh- Rod points are pixel-by— pixel measurements of iJmol in spirals, binned by the quantity on the x-axis. Green and purple points 
show tilted rings in spiral and dwarf galaxies with i?„ioi inferred from SgpR and Sjji assuming a fixed SFE {H2). We show the same data, 
binned, in Figure 17. 



In this case, we can break the SFE into two parts: star 
formation within GMCs and GMC formation, so that 



SFE = SFE (H2) ^ = SFE (H2) 



R 



■mol 



1 



(29) 



^gas Rmol 

i.e., the SFE is a product of a constant SFE (H2) and 
^moi = Sh2/Shi, which is a function of local conditions. 

We show this directly in Figure 16 and plot the same 
data, binned, in Figure 17. We plot i?,„oi — Sh2/Shi 
on the y-axis as a function of radius, S*, Ph (oc % 
and Torb- Red points show direct measurements of Rmoi 
from CO and H I assembled following the methodology 
used by Blitz & Rosolowsky (2006) to compute i?moi as 
a function of P\^: 

1. For each galaxy, we examine scatter plots to esti- 
mate a value of Ph above which our pixel-by-pixel 
measurements of I]h2 are approximately complete. 



2. Where Ph is above this limit, we measure Pmoi for 
each pixel. 

3. We sort pixels into bins based on Ph and calculate 
the average and scatter in logj^Q Pmoi for the pixels 
in each bin. 

A red point in Figure 16 corresponds to one Ph bin in 
one spiral galaxy; the x- and y-error bars indicate the 
width of the bin and the scatter in Pmoi within the bin. 
We carry out analogous procedures to compute Pmoi as 
a function of rgai, S*, and Torb- 

Because of the limited sensitivity of the CO data, these 
direct measurements of Pnioi seldom probe far below 
-Rmol = 1 and do not extend to dwarfs. Therefore we 
also use Ssfr and Shi to estimate Pmoi by assuming a 
fixed SFE (H2). For each tilted ring in both subsamples, 
we convert Ssfr into S]h2 using Equation 28. We divide 
this by the observed Shi to estimate Pmoi for that ring. 
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Fig. 17. — The data from Figure 16, binned by the quantity on the x-axis into three trends: iJniol measured pixel-by— pixel in spirals 
(red) and inferred from Egpa and Shi in (green) spiral and (purple) dwarf galaxies. Thin dashed lines show iJmol = 1 (horizontal) and our 
estimate of each quantity at the H I-to-H2 transition (vertical). Dotted lines show R^nol ^ '^fi ^ (bottom left) and R^nol oc (bottom 
right). Thick dashed lines show fits of Rniol to each quantity. 



We plot the results as green points for spirals and purple 
points for dwarf galaxies^*^. This approach — essentially 
plotting SFE (H I) in units of i?moi — allows us to es- 
timate i?moi far below the sensitivity of our CO maps. 
While this extrapolation of SFE (H2) may be aggressive, 
the quantity Ssfr/^hi must be closely related to the 
ability of H I to assemble into star-forming clouds. 

Figure 17 shows the data in Figure 16 binned by the 
quantity on the x-axis. Thin dashed lines horizontal 
show i?moi — 1, i-e., Shi = Sh2, and the value of the 
property on the x-axis that we estimate at the H I-to- 
H2 transition (§5.2 and Table 5). Dashed and dotted 
lines show fits and expectations that we discuss later in 
this section. 

In spirals, the agreement between direct measurements 

Because Ph depends on Sgas, we make a first-order correction 
to Ph in dwarf galaxies based on the estimated i?mol • 



of i?moi and estimates based on Ssfr and Shi is quite 
good. There is also general agreement between spirals 
and dwarf galaxies: the two subsamples sweep out simi- 
lar, though slightly offset, trends in all four panels. The 
magnitude of the offsets between dwarf and spiral galax- 
ies that we see in Figure 17, typically 0.2-0.3 dex, offers 
indirect evidence that differences between the subsam- 
ples — metallicity, radiation fields, spiral structure (§3.1) 
— affect cloud formation or SFE (H2) at the factor of 
~ 2-3 level. 

Figures 16 and 17 show explictly what we have already 
seen indirectly throughout §4. i?„ioi is a continuous func- 
tion of environment spanning from the H2-dominated 
(i?moi ^ 10) to H I-dominated (i?,„oi 0.1) ISM, from 
inner to outer galaxy disks, and over a wide range of 
ISM pressures. This qualitatively confirms and extends 
similar findings by Wong & Bhtz (2002) and Blitz & 
Rosolowsky (2006), which were mostly confined to the 
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inner, molecule-dominated parts of spirals. 

5.4.1. Cloud Formation Timescales 

In §2, we discuss two basic ways that i?moi might be set 
by environment. First, the timescalc to form GMCs may 
depend on local conditions. If H I and H2 are in approx- 
imate equilibrium, with the entire neutral ISM actively 
cycling between these two phases, then 

I]h2 GMC lifetime 

For constant GMC lifetimes perhaps a reasonable ex- 
tension of fixed SFE (II2) - -Rmoi is set by r (H I ^ H2). 
In §4.2.2 and §4.2.3 we saw that SFE anti-correlates with 
Tff and Torb where Shi > Sh2- If GMCs form over these 
timescales then i?moi c« % "'^ or i?nioi oc T~j!,. 

However, we found that the SFE decreased more 
steeply than one would expect if these timescales alone 
dictated iZmoh so that increasing timescale for GMC for- 
mation cannot explain all of the decline in Rmoi- Figures 
16 and 17 show this directly: dotted lines in the bottom 
two panels illustrate i?moi r^^ and i?moi « t~1^. In 
both cases the prediction is notably shallower than the 
data in both the H2 and H I-dominatcd regimes. 

5.4.2. Disk Stability Thresholds 

Of course, the entire ISM may not participate in cloud 
formation. Star formation thresholds are often invoked 
to explain the decrease in SFE between inner and outer 
galaxy disks. The amount of stable, warm H I may de- 
pend on environment, with a variable fraction of the disk 
actively cycling between H I and GMCs. This suggests a 
straightforward extension of Equation 30, 



in dwarf galaxies (see also Hunter et al. 1998a; Wong & 
Blitz 2002; Boissier et al. 2003). 

Including the effects of stellar gravity reduces stability. 
Over most of our sample, Qstars-Fgas ^ 2 with a much nar- 
rower range than Qgas (similar improvements were seen 
by Boissier et al. 2003; Yang et al. 2007). These val- 
ues arc roughly consistent with the conditions for cloud 
formation found from simulations. Li et al. (2005) find 
gas collapses where Qstars-i-gas ^1-6 and Kim & Ostriker 
(2001, 2007) find runaway instabilities where Qgas ^1-4 
(though this is Qgas and not Qstars-f-gas; for a region like 
the solar neighborhood, Kim & Ostriker 2007, argue that 
disk thickness, which tends to increase stability, approx- 
imately offsets the effect of stars on Q). 

<3stars-i-gas increases towards the central parts of spirals, 
so that although the ISM in these regions is usually dom- 
inated by H2, they appear more stable than gas near the 
H I-to-H2 transition. Hunter et al. (1998a) and Kim & 
Ostriker (2001) suggest that because of low shear, insta- 
bilities aided by magnetic fields may grow in these regions 
despite supercritical Q. Comparing to the shear thresh- 
old proposed by Himter et al. (1998a), we find some sup- 
port for this idea: at < 0.2 r25 many dwarf and spiral 
galaxies appear unstable or marginally stable. As with 
Qgas, however, Scrit,A/Sgas shows large scatter and no 
clear ability to predict the SFE. 

Thus, we find no clear evidence that disk stability at 
large scales drives the observed variations in SFE and 
i?nioi- Improved handling of second-order effects (disk 
thickness, (Tgas, -'^co, o"*, and T^) may change this pic- 
ture, but comparing our first-order analysis to expec- 
tations and simulations, disks appear marginally stable 
more or less throughout with little correlation between 
proposed thresholds and SFE. 



_ Sh2 ^, CMC lifetime 

-Kmol — ~ /TT T TT \ ^ JGMC forming j (,<jlj 

^HI T (H I ^ H2 ) 

which again balances GMC formation and destruction 
but now includes the factor /gmc forming to represent the 
fact that only a fraction of the H 1 is actively cycling 
between the molecular and atomic ISM. 

We considered three thresholds in which large-scale in- 
stabilities dictate /gmc forming Qgas, Qstars+gas, and 

shear. One would naively expect these thresholds to cor- 
respond to /gmc forming ~ 1 for supcrcritical gas and 
,/gmc forming ^ 1 for subcritical gas, yielding a step 
function in SFE or i?moi- However, we do not observe 
such relationships between thresholds and SFE (§4.3), 
which agrees with Boissier et al. (2007) who also based 
their SFR profiles on extinction-corrected FUV maps and 
found no evidence for sharp star formation cutoffs. 

If these instabilities regulate star formation but oper- 
ate below our resolution, we still expect a correspondence 
between SFE and the average threshold value, which 
should indicate what fraction of the ISM is unstable. De- 
spite this expectation, Qgas shows little correspondence 
to the SFE and almost all of our sample is stable against 
axisymmetric collapse. Kim & Ostriker (2001) and Kim 
& Ostriker (2007) discuss Qgas thresholds for the growth 
of non-axisymmetric instabilities, but these are in the 
range Qgas ^ 1-2, still lower than the typical values that 
we observe (§4.3.1). Even independent of the normaliza- 
tion, Qgas shows little relation to the SFE, particularly 



5.4.3. Cold Phase Formation 

Timescales and thresholds computed at 400 (dwarfs) 
and 800 pc (spirals) scales do not offer a simple way 
to predict i?moi- An alternative view is that physics 
on smaller scales regulates cloud formation. Comparison 
with models by Schaye (2004) suggests that a cold phase 
can form across the entire disk of most of our sample, 
which agrees with results from Wolfire et al. (2003) mod- 
eling our own Galaxy. High density, narrow-linewidth 
clouds may easily be unstable or be rendered so by the 
passage of spiral arms or supernova shocks, even where 
the ISM as a whole is subcritical. Both Schaye (2004) 
and de Blok & Walter (2006) have emphasized the effect 
of lower agas on instability and we have seen that a shift 
from the observed o-gas = 11 km s~^ to (Tgas = 3 km s~^ 
would render most gas disks in our sample iinstable or 
marginally stable (of course a proper calculation requires 
estimating the density and fraction of the mass in this 
phase as well). 

A narrow-line component is observed from high- 
velocity resolution H I observations of nearby irregular 
galaxies (Young et al. 2003; de Blok & Walter'2006), but 
an important caveat is the lack of direct evidence for 
such a component in THINGS. With ~ 2.5 or 5 km s~^ 
velocity resolution, one cannot distinguish a narrow com- 
ponent directly. Therefore, Usero et al. (2008) followed 
up on work by Braun (1997), who used the peak intensity 
along each line of sight to estimate the maximum contri- 
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TABLE 6 

FlTS^ OF R^ol = Sh2/ShI to (Ph/Po)' 



Source 


logio Po/ks 


a 




(cm- 3 K) 




spiral subsample (CO/H I) 


4.19 


0.73 


spiral subsample (SFR/H I)'' 


4.30 


0.79 


spiral subsample (combined) 


4.23 


0.80 


dwarf subsample (SFE/H 1)^ 


4.51 


1.05 


Wong & Blitz (2002) 




0.8 


Blitz & Rosolowsky (2006) 


4.54 ± 0.07 


0.92 ± 0.07 



''Over the range Rmoi = 0.1 - 10. 
'^Estimating Eh2 from Ssfr. 



bution from a cold phase and found pervasive networks 
of high brightness filaments. Usero et al. (2008) find no 

clear evidence for a cold phase traced by networks of 
high brightness filaments, suggesting that a cold phase, 
if present, is mixed with the warm phase at the THINGS 
resolution of several times ~ 100 pc. 

5.4.4. i?inoi cind Pressure 

Despite this caveat, our results offer significant circum- 
stantial support that ISM physics below our resolution 
dictates Rmoi- the lack of obvious threshold behavior, 
marginal stability of our disks, the ability of a cold phase 
to form, and the continuous variations in SFE and iJmoi 
as a function of radius, S*, and _Ph- 

In particular, the relationship between i?nioi and Ph 
has been studied before. Following theoretical work by 
Elmegrecn (1993) and Elmegreen & Parravano (1994), 
Wong & Blitz (2002) and BUtz & Rosolowsky (2006) 
showed that i^moi and Ph correlate in nearby spiral galax- 
ies (mostly at i?moi > 1) and Robertson & Kravtsov 
(2008) recently produced a similar relationship from sim- 
ulations that include cool gas and photodissociation of 
H2; they emphasize the importance of the latter to re- 
produce the observed scaling. 

The dash-dotted line in the bottom left panel of Figure 
17 shows i?moi oc P^'^, predicted by Elmegreen (1993) 
from balancing H2 formation and destruction in a model 
ISM. This is a reasonable description of dwarf galaxies, 
where wc derive a best-fit power law with index « 1.05. 
Spirals show a slightly shallower relation between Pmoi 
and Ph with best-fit power law index w 0.80. The thick 
dashed line in the bottom left panel of Figure 17 shows 
our best fit to the spiral subsample (both CO /H I and 
SFR/H I) over the range R^oi = 0.1 10. Table 6 Usts 
this fit along with fits to dwarf galaxies and the results 
of Wong & Blitz (2002) and Blitz & Rosolowsky (2006). 

The entry "spiral subsample (combined)" in Table 6 
lists the best fit power law like Equation 11 for our spiral 
subsample. This fit has an index a = 0.80 and normal- 
ization log^QPo/fce = 4.23 (this is an OLS bisector fit 
over the range 0.1 < Pmoi < 10 giving equal weight to 
each of the red and green points in Figure 16). Formally, 
the uncertainty in the fit is small because it includes a 
large number of data points. However, both logj^Q Po/fce 
and a scatter by several tcnitlis when fit to individual 
galaxies. This agrees well with a = 0.8 derived by 
Wong & Blitz (2002) and with a = 0.92 ± 0.07 obtained 
by Blitz & Rosolowsky (2006) given the uncertainties. 
Fitting the dwarf subsample in the same manner yields 
logio Po/^B = 4.51, the pressure at the H I-to-H2 transi- 



tion. This is 0.2-0.3 dex higher than \ogiQ Po/fce = 4.23 
in spirals, suggesting that at the same pressure (density) 
GMC/H2 formation in our dwarf subsample is a factor 
of 2 less efficient than in spirals. 

5.4.5. Pmoi and Environment 

The fits between Pmoi and Ph in Table 6 are reasonable 
descriptions of the data, but do not represent a "smoking 
gun" regarding the underlying physics; radius, S,, Ph, 
and Torb are all covariant and each could be used to pre- 
dict Pmoi with reasonable accuracy in spirals. Therefore 
we close our discussion by noting a set of four scaling re- 
lations between Pmoi and environment that describe our 
spiral subsample 

Rmoi = 10.6 exp (-rgai/0.21 r25) (32) 
Pmoi = S*/81 M0 pc-2 (33) 

Pmol = (Ph/1.7 X 10^ cm-3 K fcB)° ' (34) 

Pmol = (Torb/1-8 X 10^ yr) "^ ^ (35) 

these appear as thick dashed lines in Figure 17. 

In particular, we stress the relationship between Pmoi 
and (see also Figure 3). This has several possible in- 
terpretations, the most simple of which is that stars form 
where they have formed in the past. There are physical 
reasons to think relationship may be causal, however. 
Considering a similar finding in dwarf irregular galax- 
ies. Hunter et al. (1998a) suggestcxl that stellar feedback 
may play a critical role in triggering cloud formation. 
Recently the importance of the stellar potential well has 
been highlighted, either to triggering large-scale instabil- 
ities (Li et al. 2005, 2006; Yang et al. 2007) or in bringing 
gas to high densities in order for small-scale physics to 
operate more effectively (Elmegreen 1993; Elmegreen & 
Parravano 1994; Wong & Blitz 2002; Blitz & Rosolowsky 
2004, 2006). 

5.5. A Note on Systematica: Xqo, <^ga.s, c*; 

In this paper, we work "to first order," using the sim- 
plest well-motivated assumptions to convert observations 
to physical quantities. These assumptions are described 
in §3 and Appendices A - D. These are not always unique 
and here we note differences with the literature and the 
effect that they may have on our analysis. 

Xco- In spirals, we adopt a fixed Xqo = 2 x 
10^0 cm-2 (Kkms-i)-i. Wong & Blitz (2002) and 
Blitz & Rosolowsky (2006) adopt the same value. Ken- 
nicutt (1989), Kennicutt (1998a), Martin & Kenni- 
cutt (2001), and Kennicutt et al. (2007) also use a 
fixed value, but a slightly higher one, Xco = 2.8 x 
10^0 cm-2 (K km s''^)-'^. Boissier et al. (2003) test the 
effects of a metallicity-dependent Xco that tends to 
yield lower Sh2 than our values in the inner parts of 
spirals, but higher in the outer parts. 

Variations in the normalization of Xco will affect the 
location of the H I to H2 transition and the value of 
SFE (H2), but not the observations of fixed SFE (H2) 
or steadily varying Pmoi- A strong dependence of Xco 
on environment in spirals would affect many of our re- 
sults, but leave the basic observation of environment- 
dependent SFE (H I) intact. Variations in the CO J = 
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2^1/1—^0 line ratio (Appendix A) will manifest as 
changes in Xco- 

Cgas-' We adopt cigas = 11 km. based on the 
THINGS second moment maps (Appendix B). This is 
almost twice the commonly used Ugas = 6 km s~^ (Ken- 
nicutt 1989, 1998a; Martin & Kennicutt 2001; Boissier et 
al. 2003) and also higher than cTgas ~ 8 km s~^, adopted 
by Blitz & Rosolowsky (2006). Wc emphasize the impor- 
tance of CTgas to the stability analysis in §4.3.1; observa- 
tions with velocity and spatial resolution capable of dis- 
tenangling different H I components and a multi-phase 
analysis arc needed to move forward on this topic. 

CTh.; We assume an isothermal stellar disk with a fixed 
scale height, as do Boissier et al. (2003) and Blitz & 
Rosolowsky (2006). Wong & Blitz (2002) and Yang 
et al. (2007) assume a fixed stellar velocity dispersion. 
This has a moderate effect on Ph and a strong effect on 

Qstars-)-gas (§4-3.2). 

Tf .• We adopt Tf = 0.5Mq/Lq^k (Appendix C), 
consistent with our adopted IMF and Bell et al. (2003), 
the same value used by Blitz & Rosolowsky (2004, 2006). 
Tf directly affects S,, Ph, re, and Qstars-i-gas- It may 
vary by ~ 30% both within and among galaxies (Bell 
& de Jong 2001), with larger variations in the bluest 
galaxies or from changes to the assumed IMF (Bell et al. 
2003). 

Star formation rate tracer: We use FUV+24/xm to es- 
timate recent Esfr, discussed in detail in Appendix D. 
This is similar to Boissier et al. (2007) but in contrast 

to Kennicutt (1989), Kennicutt (1998a), Martin & Ken- 
nicutt (2001), Wong & Blitz (2002), and Boissier et al. 
(2003), who each used Ha emission with various extinc- 
tion corrections. Boissier et al. (2007) considered differ- 
ences between between Ha and FUV profiles in detail, 
suggesting that stochasticity leads Ha to show signs of 
knees and turnoffs while FUV remains smooth. We work 
pixel by pixel and in radial profile (similar to Martin & 
Kennicutt 2001; Wong & Blitz 2002; Boissier et al. 2003; 
Blitz & Rosolowsky 2006) rather than attempting to iso- 
late individual star forming regions (e.g., Kennicutt et al. 
2007). Both differences mean that we measure "recent" 
rather than "present" Ssfr, which may account for some 
of the smoothness in the trends seen in §4. 

6. STAR FORMATION RECIPES 

As a final exercise, we compare our galaxies to simple 
star formation recipes based on the laws and thresholds 
discussed in §2 and normalized to the H I-to-H2 transi- 
tion in spirals (§5.2). We predict the SFE in this way: 

1. We assume SFE (H2)= 5.25 x 10-^° yr"! 

2. We calculate Rmoi = ^H2/^m either: 

(a) By setting Pmoi = rsfi/rs or Torb,o/Torb- 

(b) From our fits of i?moi to radius, E*, Ph, and 
Torb in spirals galaxies (Equation 32). 

3. We derive Ssfr from SFE (H2), Pmoi, and Ehi- 

4. We calculate the predicted SFE, dividing Ssfr by 
Sgas in spirals and Shi in dwarf galaxies. 

5. We combine Rmoi t^^ or r^j.^ with thresholds. 
In subcritical areas, we set SFE = 5 x 10~^^ yr~^, 



roughly the observed value at r25 in both subsam- 
ples. 

Figure 18 illustrates the procedure for Pmoi oc r^^ and 
the Qstars-i-gas threshold in the spiral galaxy NGC 3184. 

We set Ts,o and Torb,o equal to the timescale at the 
H I-to-H2 transition in spirals (§5.2, Table 5), i.e., we 
predict Pmoi using the dotted lines in Figure 16. The pre- 
dictions will therefore intersect our data where Pmoi = 
Eh2/Shi » 1 in spirals. 

We adopt the same approach to normalize thresholds. 
For shear and Qstars+gas, we define the boundary be- 
tween supercritical and subcritical data as 2.3 and 1.6, 
respectively, approximately the values at the H I to H2 
transition in spirals. For the Schaye (2004) cold phase 
threshold we use a critical value of 1. 

We implement thresholds pixel-by-pixel and present 
our results in radial average. Within a tilted ring, some 
lines of sight can be supercritical and some can be sub- 
critical, allowing the threshold to damp the average SFE 
in a ring without setting it to the minimum value. 

The choice to normalize the recipes for both dwarf and 
spiral galaxies using values measured for spirals is meant 
to highlight differences between the subsamples. 

6.1. Results 

Figure 19 shows the results of these calculations. The 
observed SFE as a function of radius appears as a shaded 
gray region (based on Figure 1). Radial profiles of SFE 
compiled from predictions appear in color (these fol- 
low the same methodology used to make the bins in 
Figure 1). The top row shows results for spiral (lejf) 
and dwarf {right) galaxies setting Pmoi oc combined 
with several thresholds, the middle row shows results for 
-Rmoi oc combined with the same stable of thresholds, 
and the bottom row shows Pmoi set by fits to radius, S*, 
Ph, and Torb- 

Figure 19 illustrates much of what we saw in §4 and 5. 
First, adopting fixed SFE (H2) ensures that we match the 
observed SFE with reasonable accuracy in the inner parts 
of spirals regardless of how we predict Pmoi- Using fits to 
predict Pmoi (bottom row) offers a small refinement over 
the timescales in this regime, but as long as Pmoi ^ 1 
then SFE ~ SFE (H2)- As a result, the available gas 
reservoir sets the SFR in this regime. 

Setting Pmoi % ^ or r^j'^ smoothly damps the SFE 
with increasing radius, but not by enough to match ob- 
servations. Without a threshold, Tff and Torb overpredict 
the SFE at large radii in spiral galaxies and at r > 0.4 r25 
in dwarf galaxies (black bins, top two rows). 

Thresholds damp the SFE with mixed success (green, 
magenta, and blue bins in the top two panels). Each 
somewhat lowers the SFE in the outer parts of galax- 
ies. In the process, however, both Qstars+gas and shear 
predict suppressed star formation at low or intermediate 
radii in both dwarf and spiral galaxies, areas where we 
observe ongoing star formation (the vertical error bars 
show that the 50% range includes completely subcritical 
galaxies in both cases). Wc saw in §4.3 that the radial 
variation in these thresholds is often less than the scat- 
ter among galaxies at a given radius and that the step 
function behavior that we implement here is not clear in 
our data. 

The Schaye (2004) threshold predicts that a cold phase 
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Fig. 18. — How we predict Sgppi, illustrated for NGC 3184. VVc calculate the threshold value {left, here Qstars+gas)i identify supercritical 
areas (solid contour). In parallel, we estimate Rmol = Sh2/Shi (middle), here from Rj^ol "^g • We combine these with the assumption 
of a fixed SFE {H2) and a (fixed) low SEE in subcritical areas to predict Ssfr (right). 



can form almost everywhere in our sample and so only 
comes into play in the outer parts of spirals and in dwarf 
galaxies, where it damps the predicted SFE, but not by 
enough to match observations. 

The bottom left panel shows that the fits mostly do 
a good job of reproducing the SFE in spirals, which is 
expected because they arc fits to these data. 

The same fits (to spirals) yield mixed results when ap- 
plied to dwarf galaxies. The fits to and Torb show 
very large scatter and fits to radius, Ph, and Torb all 
overpredict SFE by varying amounts (similar discrepan- 
cies are evident comparing spirals and dwarfs in the top 
two panels). The scaling relations relating Rmoi to envi- 
ronment in spirals apparently do not apply perfectly to 
dwarfs. Likely drivers for the discrepancy are the lower 
abundance of metals and dust and more intense radiation 
fields, which afi'ect phase balance in the ISM and the rate 
of H2 formation and destruction. Focusing on the pres- 
sure fit (green), we can phrase the observation this way: 
for the same pressure (density), cloud formation in our 
dwarf subsample is suppressed relative to that in spirals 
by a factor of ~ 2. 

7. SUMMARY 

We combine THINGS, SINGS, the GALEX NGS, 
HERACLES, and BIMA SONG to study what sets the 
star formation efiiciency in 12 nearby spirals and 11 
nearby dwarf galaxies. 

We use these data to estimate the star formation rate 
surface density, gas kinematics, and the mass surface 
densities of HI, H2, and stars (Appendices A and C). 
To trace recent star formation, we use a linear combina- 
tion of GALEX FUV and Spitzer 24/im (Appendix D). 
We suggest that this combination represent a useful tool 
given the outstanding legacy data sets now available from 
these two observatories (e.g., SINGS and the GALEX 
NGS). 

We focus on the star formation efiiciency (SFE), 

SsFR/Sgas, and the H2-to-H I ratio, i?moi- These 
quantities remove the basic scaling between gas and 
SFR, allowing us to focus on whore gas forms stars 
quickly /efficiently (SFE) and the phase of the neutral 
ISM (i?moi)- We measure the SFE out to ~ 1.2 r25, com- 
pare it to a series of variables posited to influence star 
formation, and test the ability of several predictions to 
reproduce the observed SFE. 



7.1. Structure of Our Typical Spiral and Dwarf Galaxy 

We deliberately avoid discussing individual galaxies in 
the main text (these data appear in Appendices E and 
F). Instead, we study "stacked" versions of a spiral and 
dwarf galaxy. We sketch their basic structure here. 

The spiral galaxy has a roughly constant distribution of 
H I, Shi ^ 6 M0 pc^^ out to ^ r25. H I surface densities 
seldom exceed Shi ~ 10 Mq pc~^; gas in excess of this 
surface density tends to be molecular. Wo observe no 
analogous saturation in Eh2, finding Sh2 ^ 100 Mq pc~^ 
in the very central parts of many galaxies. 

Molecular gas, star formation, and stellar surface 
density all decline with nearly equal exponential scale 
lengths, ^0.2 r25, giving the appearance of a long lived 
star-forming disk embedded in a sea of H I. The ISM is 
mostly H2 within ~ 0.5 r25 and where S* > 80 Mq. 

Over a wide range of conditions the SFR per unit H2, 
SFE (H2), = 5.25 ± 2.5 x 10^1° yr^^ at scales of 800 pc. 
This is a "limiting efficiency" in the sense that we do not 
observe the average SFE in spirals to climb above this 
value. Where the ISM is mostly H I, the SFE is lower 
than this limiting value and declines radially with an ex- 
ponential scale length ~ 0.2-0.25 r25. In this regime, the 
star formation rate per unit stellar mass remains nearly 
fixed at a value about twice the cosmologically average 
rate (i.e., the stellar assembly time is ~ twice the Hubble 
time). 

Dwarf galaxies also exhibit fiat H I distributions, de- 
clining SFE with increasing radius, and a nearly constant 
stellar assembly time. Normalized to r25, the scale length 
of the decline in the SFE is identical to that observed 
in spirals within the uncertainties. The stellar assem- 
bly time is half that found in spirals, corresponding to 
roughly a Hubble time. Dwarfs exhibit only the crudest 
relationship between Esfr and Ehi and, as a result, 
is a much better predictor of the SFR than Sgas (in good 
agreement with Hunter et al. 1998a). The lack of a clear 
relationship between Ssfr and Sgas is at least partially 
due to an incomplete census of the ISM: conditions in the 
central parts of dwarf galaxies often match those where 
we find H2 in spirals and in these same regions the SFE 
is (unexpectedly) higher than we observe anywhere in 
spirals (where H2 is included). 

7.2. Conclusions for Specific Laws and Thresholds 



30 



Leroy et al. 




Observations ■ 

Threshold: 
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shear 4- 
cold phuae + 




0.2 0.4 0.6 O.B 1.0 
Celaclocenlric Radius [fsj] 



Observations ■' 
Threshold: 
none + 
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cold phase + 




0.2 0.4 0.6 O.B 1.0 
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0,0 O.Z 0.4 0.6 O.S 1.0 I,Z 0,0 0.2 0.4 0.6 O.B 1.0 1.2 

Galactooentric Radius [rjs] Galactoeentric Radius [rsj] 



Fig. 19. — Comparison of predicted (color bins) to observed (gray region) SFE in spiral (left) and dwarf (right) galaxies. We adopt fixed 
SFE (H2) and predict i^mol from (top row) and (middle panel) combined with thresholds. We also show four fits of Rmol to other 
quantities in spirals {bottom row). The dotted horizontal line in the top two rows shows the SFE that we adopt for subcritical data. 



The Star Formation Efficiency in Nearby Galaxies 



31 



Wc compare the observed SFE to proposed star for- 
mation laws and thresholds described in §2. For star 
formation laws we find: 

• The SFE varies dramatically over a small range of 
Ehi and very little with changing Sh2- Therefore, 

the disk free-fall time for a fixed scale height disk 
or any other weak dependence of SFE on Sgas is of 
little use to predict the SFE (§4.2.1). 

• The disk free-fall time accounting for a changing 
scale height, Tfj , correlates with both SFE and iJmoi 
(§4.2.2). Sotting SFE proportional to Tff broadly 
captures the drop in SFE in spirals, but predicts 
variations in SFE (H2 ) that we do not observe and 
is a poor match to dwarf galaxies. Taking rg- to be 
the relevant timescale for H I to form GMCs (i.e., 
-Rmoi oc Tff ^), fails to capture the full drop in the 
SFE in either subsample. 

• The orbital timescale, Torb, also correlates with 
both SFE and i?moi, but in outer spirals and dwarf 
galaxies both SFE and Rmo\ drop faster than Torb 
increases (§4.2.3). As with rs, Torb alone cannot 
describe cloud or star formation in our sample. 

• In spirals, we observe no clear relationship between 
SFE (H2) and the logarithmic derivative of the ro- 
tation curve, (5 (§4.2.3). In dwarf galaxies, SFE 
correlates with /3. Both observations are contrary 
to the anti-correlation between SFE and /3 expected 
if cloud-cloud collisions set the SFE (Tan 2000). 

• Fixed GMC efficiency appears to be a good de- 
scription of our spiral subsample (§4.1 and Bigicl 
et al. 2008). SFE (H2) is constant as a function of 
a range of environmental parameters. This obser- 
vation applies only to the disks of spiral galaxies, 
not starbursts or low metallicity dwarf galaxies. 

• Wc observe a correspondence between hydrostatic 
pressure and ISM phase (§4.2.2 and §5.4.4). In spi- 
rals our results are consistent with previous work 
(Wong & Blitz 2002; Blitz & Rosolowsky 2006). In 
dwarf galaxies and the outer parts of spirals, in- 
ferring i?nioi from SFE (H I) yields results roughly 
consistent with predictions by Elmegreen (1993). 

For thresholds we find: 

• Despite a suggestion of increased stability at large 
radii in spirals, there is no clear relation between 
Qgas — which measures stability against axisym- 
metric collapse due to self-gravity in the gas disk 
alone — and SFE. Most regions are quite stable 
and Qgas has large scatter, even appearing weakly 
anii-correlated with the SFE in dwarfs (§4.3.1). 

• When the cfi^ects of stars are included, most disks 
are only marginally stable: Qstars+gas (Rafikov 
2001), which measures gravitational instability in 
a disk of gas and stars, lies mostly in the nar- 
row range 1.3-2.5, increasing slightly towards the 
centers and edges of galaxies. We emphasize that 
adopted parameters — ^CO) Cgas, T^, and a^, — 
strongly affect both Qgas and Qstars-Fgas (§4.3.2). 



• The ability of instabilities to survive competition 
with shear (Himter et al. 1998a) shows the same 
large scatter and high stability as Qgas in the outer 
disks of spirals, but identifies most areas in dwarf 
galaxies and inner spirals as only marginally stable, 
an improvement over Qgas (§4.3.3). 

• Most areas in both dwarf and spiral galaxies meets 
the condition needed for a cold phase to form 
(§4.3.4) (Schaye 2004). Regions that do not meet 
this criterion tend to come from outer disks and 
have low SFE. Because this criterion is met over 
such a large area, it is of little use on its own to 
predict variation in the SFE within galaxy disks. 

Finally, we distinguish three different critical surface 
densities. First, in spirals Egas ~ 14 M0 pc~^ at the 
H I-to-H2 transition. We find no evidence that this is 
a real threshold for cloud formation: i?,noi = Sh2/Shi 
varies continuously across i?inoi = 1 as a function of other 
quantities. However, it is useful to predict the SFE, 
which will be nearly constant above this Sgas- A related 
(but not identical) value, Shi ~ 10 M0 pc~^, is the sur- 
face density at which H I "saturates." Gas in excess of 
this surface density is in the molecular phase (Martin & 
Kennicutt 2001; Wong & Blitz 2002; Bigiel et al. 2008). 
This presumably drives the; observation that most vigor- 
ous star formation takes place where Shi ^ 10 Mq pc~^ 
(e.g., Skillman 1987). Last, lower values, S 

gas ^~ 

4 Mq pc~2 (e.g. Kennicutt 1989; Schaye 2004), may cor- 
respond to the edge of the star-forming disk. At our 

resolution such values are relatively rare inside 1.2 r25 
and we draw no conclusion regarding whether this "outer 
disk threshold" corresponds to a real shift in the mode 
of star formation. 

7.3. General Conclusions 
Our general conclusions are: 

1. In the disks of spiral galaxies, the SFE of II2 is 
roughly constant as a function of: galactocentric 
radius, S,, Sgas, ^h, Torb, Qgas, and P (§5.1). 
This fixed SFE (H2)= 5.25 ± 2.5 x 10"!° yr"! 
(tdcp(H2) = 1.9 X 10^ yr) sets the SFE of total 
gas across the II2 dominated inner parts (rgai < 
0.5 r25) of spiral galaxies. 

2. In spiral galaxies, the transition between a mostly- 
H I and a mostly-H2 ISM is a well-defined func- 
tion of local conditions (§5.2). It occurs at a 
characteristic radius (0.43 ± 0.18 r25), S, (81 ± 
25 Mq pc-2), Sgas (14 ± 6 Mq pc-^), Ph (2.3 ± 
1.5 X W fee cm--"' K), and Torb (1.8±0.4 x 10« yr). 

3. We find indirect evidence for abundant H2 in 
the central parts of many dwarf galaxies, where 

SFE (H I) exceeds SFE (H2) found in spirals. The 
simplest explanation is that H2 accounts for a 
significant fraction of the ISM along these lines 
of sight (§4.1.1 and §5.3). The implied central 
Sh2/Shi is ~ 2.5 with Sh2 = Shi at ~ 0.25 r25- 

4. Where Shi > Sh2 — in the outer parts of spirals 
and throughout dwarf galaxies (by assumption) — 
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we observe the SFE to decline steadily with increas- 
ing radius, with scale length ^ 0.2 0.25 r25 in both 
subsamples (§4.1). We also observe a decline in 
SFE with decreasing S*, decreasing Ph, and in- 
creasing Torb, which are all covariant with radius. 

5. Where Shi > Sh2, we find little relation between 
SFE and Sgas (§4.1.3) but a strong relationship be- 
tween SFE and (§4.1.2). The simplest explana- 
tion is that present day star formation roughly fol- 
lows past star formation. A more aggressive inter- 
pretation is that the stellar potential well or feed- 
back are critical to bring gas to high densities. 

6. The H2 to H I ratio, i?inoi = Sh2/Shi, and by ex- 
tension cloud formation, depends strongly on envi- 
ronment. Rmoi correlates with radius, Ph, ts, Torb; 
and S* in spirals. We find corresponding corre- 
lations between these quantities and Esfr/Shi, a 
proxy for the efficiency of cloud formation in dwarfs 
and the outer parts of spirals. At our resolution, 
-Rmoi appears to be a continuous function of envi- 
ronment from the H I-dominated (i?moi ^ 0.1) to 
H2-dominated (Pmoi ~ 10) regime (§5.4). 

7. The variation in P,noi is too strong to be repro- 
duced only by varying Torb or Tg (§4.2 and 5.4.1). 
Physics other than these timescalcs must also play 
an important role in cloud formation (points 8 - 
11). 

8. Thresholds for large scale stability do not offer an 
obvious way to predict Pmoi- We find no clear re- 
lationship (continuous or step-function) between 
SFE and Qgas, Qstars+gas, Or the shear threshold. 
The threshold values we find suggest disks that are 
stable or marginally stable throughout once the ef- 
fects of stars are included (§4.3 and 5.4.2). 

9. We derive a power law relationship between i?inoi 
and hydrostatic pressure (Elmegrcen 1989) that is 
roughly consistent with expectations by Elmegreen 
(1993), observations by Wong & Blitz (2002) and 
Blitz & Rosolowsky (2006), and simulations by 
Robertson & Kravtsov (2008). In its simplest form, 
this is a variation on the classical Schmidt law, i.e., 
-Rmoi set by gas volume density (§4.2 and 5.4.4). 



10. Power law fits of i?,noi to Ph (rff), radius, Torbj 
and reproduce observed SFE reasonably in spi- 
ral galaxies but yield large scatter or higher-than- 
expected SFE in the outer parts of dwarf galax- 
ies, offering indirect evidence that the differences 
between our two subsamples — metallicity (dust), 
radiation field, and strong spiral shocks — play a 
role in setting these relations (§5.4.5 and 6). 



11. Our data do not identify a unique driver for the 
SFE, but suggest that ISM physics below our res- 
olution — balance between warm and cold H I 
phases, H2 formation, and perhaps shocks and tur- 
bulent fluctuations driven by stellar feedback — 
govern the ability of the ISM to form GMCs out of 
marginally stable galaxy disks (§5.4.5). 
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APPENDIX 

These appendices describe how we assemble the database of radial profiles and maps that are used in the main text. 
We discuss the data and methods that wc use to derive gas surface densities (Appendix A), kinematics (Appendix B), 
stellar surface densities (Appendix C), and star formation rate surface densities (Appendix D). Finally, we present a 
table containing radial profiles of key quantities (Appendix E) and an atlas showing maps, profiles, and basic results 
for each galaxy (Appendix F) . 

MAPS OF H I AND H2 SURFACE DENSITY 

Shi from THINGS 21cm Maps 

The H I Nearby Galaxy Survey (THINGS, Walter et al. 2008) mapped 21-cm line emission from all of our sample 
galaxies using the Very Large Array. We calculate atomic gas mass surface density, Shi, from natural- weighted data 
that have mean angular resolution 11" and mean velocity resolution 5 km s~^. THINGS includes data from the most 
compact VLA configuration and therefore comfortably recovers extended structure (up to 15') in our sources. At 
30" resolution, THINGS maps are sensitive to Shi low as ^ 0.5 Mq pc^^; here wc adopt a working sensitivity of 
Shi = 1 pc~^. In practice the sensitivity and field of view of the THINGS maps are sufficient to measure Shi to 
> r25 in almost every galaxy. For detailed description and presentation of THINGS, we refer the reader to Walter et 
al. (2008). 

To convert from integrated intensity to Shi we use 

Shi [Mq pc"2] = 0.020 cos i /21cm [K km s~^] . (Al) 
which accounts for inclination and includes a factor of 1.36 to reflect the presence of helium. 

Sh2 from HERACLES (IRAM 30-m) and BIMA SONG CO Maps 

We estimate the surface density of molecular hydrogen, Sh2, from CO emission, the most commonly used tracer of 
H2. Along with Bigiel et al. (2008), this study presents the first scientific results from HERACLES, a large project 
that used the HERA focal plane array (Schuster et al. 2004) on the IRAM 30-m telescope to map CO J = 2 ^ 1 
emission from the full optical disk in 18 THINGS galaxies (Leroy et al. 2008). These data have an angular resolution 
of 11" and a velocity resolution of 2.6 km s~^. The typical noise in an individual channel map is 40-80 mK, yielding 
(masked) integrated intensity maps that are sensitive to Sh2 ^ 4 Mq pc~^ at our working resolution and adopted 
conversion factor. 

HERA maps are not available for NGC 3627 and NGC 5194. In these galaxies, we use CO J = 1 ^ maps from 
the BIMA Survey of Nearby Galaxies (BIMA SONG, Heifer et al. 2003) to estimate Sh2- These data have angular 
resolution 7" and include zero-spacing data from the Kitt Peak 12m, ensuring sensitivity to extended structure. 

We derive Sh2 from integrated CO intensity, Iqo by adopting a constant CO-to-H2 conversion factor, Xqo = 
2 X 10^*^ cm^^ (K km s^^)^^. Based on comparison to 7-ray and FIR observations, this value is appropriate in the 
Solar Neighborhood (Strong & Mattox 1996; Dame et al. 2001). For CO J = 1 ^ emission, the conversion to Sh2 is 

Sh2 [Mq pc-2] = 4.4 cosi /co (1^0) [K km s"^] . (A2) 

To relate CO J = 2 ^ 1 to CO J = 1 ^ intensity, we further assume a line ratio of /co(2 ^ 1) = 0.8 /co(l ^ 0). 
Based on direct comparison of HERACLES and previous surveys, thisis a typical value in our sample (Leroy et al. 

2008) and is intermediate in the range 0.6 1.0) observed for the Milky Way and other spiral galaxies (, e.g., 

Braine et al. 1993; Sawada et al. 2001; Schuster et al. 2007). Thus, for the HERACLES maps we derive Sh2 via 

Sh2 [M0 pc-2] = 5.5 cosi Jco (2^1) [K km s"^] . (A3) 

The CO-to-}l2 Conversion Factor 

The CO-to-H2 conversion factor is presumably a source of significant systematic uncertainty in Sh2 • -^co almost 
certainly varies: it is likely to be lower than Galactic (yielding lower Sh2) in overwhelmingly molecular, heavily excited 
regions; it is likely to be higher (yielding higher Sh2) in regions with low dust content and intense radiation fields, 
such as dwarf irregular galaxies. There is compelling evidence for both senses of variation, but it is our assessment 
that no reliable calibration of Xco as a function of metallicity, radiation field, and Sh2 yet exists. A useful calibration 
must reflect all of these quantities, which all affect Xco and are not universally covariant. 

In light of this uncertainty, our approach is: 1) to treat Xco as unknown in low mass, low-metallicity galaxies, 
where different approaches to measure Sh2 yield results that differ by an order of magnitude or more and 2) to assume 
that variations in Xco within spiral galaxies are relatively small. The second point might be expected based on 
theoretical modeling of GMCs (Wolfire et al. 1993) and the observed uniformity of GMC properties across a wide 
range of environments (Bolatto et al. 2008). We emphasize that even if present, the most extreme variations are likely 
to contribute primarily to the central resolution element, which is not the focus of this study, and the far outer disk, 
where Sh2 is not the dominant mass component. 

Variations aside, estimates of "typical" values of Xco in the Milky Way and other spiral galaxies span the range 

1.5-4 X 10^° cm~2 (K km s~^)~-^ (e.g. Blitz et al. 2007; Draine et al. 2007, in addition to the references already 
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Fig. B1. — lUiostration of our rotation curve treatment. Light gray profiles show median rotation velocity and scatter measured directly 
from the THINGS first moment maps; dark gray profiles and scatter show the higher quality (de Blok et al. 2008) rotation curves. The thick 
black lines show the fit that we use to approximate the rotation curve. This simple function (Equation Bl) does a good job of capturing 
both the steadily rising rotation curves typical of dwarf galaxies (e.g., IC 2574, left panel) and the rapidly rising then flat curves seen in 
more massive spiral galaxies (e.g., NGC 3198, right panel). 

given). The choice of mean Xqo within this range can have a large impact on, e.g., assessing gravitational stability or 
conditions at the H2-to-H I transition. We refer the reader to Boissier et al. (2003) for a quantitative exploration of 
how different assumptions regarding Xqo affect a stability analysis. 

Masking the H I and CO Data Cubes 

The H I and CO data cubes have large bandwidth, only a small part of which contains signal of the spectral line. 
In order to produce integrated intensity maps with good signal-to-noisc ratio, we blank signal-free regions of the 
H I and CO cubes. Walter et al. (2008) describe this process for THINGS. We apply an analogous procedure to 
the HERACLES and BIMA SONG data. We convolve the cubes to 30" resolution, identify regions with significant 
emission, and then blank the original data cubes outside these regions. We integrate these blanked cubes to create 
intensity maps. For HERACLES, we require Iqq > 2(Trms in 3 consecutive (2.6 km s^^) channels at 30" resolution. 
Note that our use of masking drives the small 10%) numerical differences with the HERACLES survey paper (which 
uses a different approach to create integrated intensity maps). For BIMA SONG, we require either Iqo > 3ijrms in a 
single (10 km s~^) channel at 30" resolution or Ico > 20-rms in consecutive velocity channels, similar to the original 
masking by Heifer et al. (2003). In both cases we consider only CO emission within <^ 100 km s~^ of the mean H I 
velocity. 

KINEMATICS 

Rotation Curves from THINGS 
We approximate all galaxies to have rotation curves with the following functional form (Boissier et al. 2003): 



1 - exp ( - — 

tfiat 



(Bl) 



where v^ot is the circular rotation speed of the galaxy at a radius r and Ufiat and /fiat arc free parameters that represent 
the velocity at which the rotation curve is flat and the length scale over which it approaches this velocity. For a 
continuously rising rotation curve, common for low-mass galaxies, we expect large /flat, while the almost flat rotation 
curves of massive spiral galaxies will have small /flat and then remain nearly constant at Uflat . 

In most cases. Equation Bl captures the basic behavior of the rotation curve well. Small scale variations are lost, 
but these may be due to streaming motions near spiral arms or warps in the gas disk as easily as real variations in the 
circular velocdty. On the other hand. Equation Bl offers the distinct advantage of having a smooth, analytic derivative. 
Our analysis uses the rotation curve to estimate the orbital timescale, shear, and coriolis force (see §2). The former is 
quite reasonably captured by Equation Bl and the latter two depend critically on the derivative of the rotation curve 

P=d\0g U(rgal)/c/l0g Tgal. 

For each galaxy, we derive Vflat and /flat from a non- linear least squares fit using Equation Bl and profiles of v^ot 
measured from the THINGS data cubes. We calculate v^ot from the intensity-weighted first moment, Vr, via 
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Vr - V, 



'sys 



(B2) 



f rot = 



e 



Sim cos 



Here Vgys is the systemic velocity, i is the inclination, and 6 is the azimuthal angle relative to the receding major axis 
measured in the plane of the galaxy. Wc calculate maps of Uj-ot and then convert these into profiles of the median and 
la scatter in v^ot within 60° of the major axis in a series 5"-wide tilted rings. We fit Equation Bl to the profile of 
median Vrot weighted by the scatter in that ring. 

For many of our galaxies, high quality rotation curves are available from the analysis of de Blok et al. (2008, see 
Table 2). Wherever possible, we include these in our fit with very high weight, so that tlioy drive the best-fit Vfiat and 
^flat for these galaxies. For the 7 low inclination galaxies in our sample that are not part of the study by de Blok et 
al. (2008) (see Column 4 of Table 2) we only fit the first moment data. 

Figure Bl shows examples of this procedure for two galaxies: the dwarf irregular IC 2574, which has a steadily rising 
rotation curve, and the spiral galaxy NGC 3198, which has a quickly rising rotation curve that remains flat over most 
of the disk. We plot Vrot and associated scatter, the de Blok et al. (2008) rotation curve, and the best-fit version 
of Equation Bl. The best -fit values of fflat and Zflat for all galaxies arc given in Table 4; for the three galaxies that 
overlap the sample of Boissier et al. (2003), we match their fitted parameters well. 

The dynamics of the irregular galaxies NGC 3077 and NGC 4449 are not well-described by Equation Bl; the former 
is disturbed by an ongoing interaction with M81 and the latter has a counter-rotating core, perhaps due to a recent 
interaction (Hunter et al. 1998b). We neglect both galaxies in the kinematic analyses. 



Throughout this paper, we adopt a single gas velocity dispersion, o-gas = 11 km s"^. This is typical of the outer 

(H I-dominatcd) parts of THINGS galaxies and agrees well with values derived by Tamburro et al. (in prep.), who are 
conducting a thorough study of Cgas in THINGS. The left panel in Figure CI motivates this choice. Wc plot the median 
and 1(7 range of CTgas over the range 0.5 r25 1.0 for each galaxy in THINGS as a function of the inclination of the 
galaxy. We restrict ourselves to the outer disk because over this regime H I usually dominates the ISM. This figure 
shows that a fixed o-gas = 11 ± 3 km s~^ is a good description of the outer disk for galaxies with i < 60°; variations 
both within and among galaxies are comparatively small, typically 25%. On the other hand, highly inclined galaxies 
show large scatter and systematically high velocity dispersions, likely because the velocity dispersion is significantly 
affected by projection effects. 

Variations in the gas velocity dispersion inside 0.5 r25 could be expected to take two forms: cTgas in the warm neutral 
medium may increase in regions of active star formation due to stellar feedback (e.g. Dib et al. 2006) and the fraction 
of gas in a narrow-line width (cold) H I phase may increase towards the centers of galaxies (e.g. Schaye 2004). The first 
effect may be observed in THINGS: the second moment maps show a gradual increase in Cgas from the outskirts to the 
centers of galaxies. The second effect can, in principle, be observed using 21-cm line observations (de Blok & Walter 
2006), but doing so is very challenging, requiring better spatial and velocity resolution and a higher signal-to-noise 
ratio than is achieved in most THINGS targets. Further, we know that a large fraction of the ISM is H2 in the central 
parts of our spiral galaxies, making it even more complicated to interpret measurements based only on H 1. Because 
measuring the detailed behavior of o-gas inside ~ 0.5 r25 is beyond the limit of our current data, and because Ogas 
varies only gradually in the outer parts of galaxies, we adopt a fixed Ogas (an almost universal approach in this field, 
following Kennicutt 1989; Hunter et al. 1998a; Martin & Kennicutt 2001; Wong & Blitz 2002; Boissier et al. 2003). 



Direct measurements of the stellar velocity dispersion, o, , across the disks of nearby galaxies are extremely scarce. 
In lieu of such observations for our sample, we make four assumptions to estimate o*. First, we assume that the 
exponential stellar scale height, /i*, of a galaxy does not vary with radius. This is generally observed for edge-on 
disk galaxies (van der Kruit & Searle 1981; de Grijs & Peletier 1997; Kregel et al. 2002). Second, we assume that /i* 
is related to the stellar scale length, by U/h^ = 7.3 ± 2.2, the average flattening ratio measured by Kregel et al. 
(2002). Because we measure this yields an estimate of /i*. Third, we assume that our disks are isothermal in the 

2;-direction, so that hydrostatic equilibrium yields /i* = 1/2 (oj _j/27rGp*)°'^ (van der Kruit & Searle 1981), where 
p* is the midplane stellar volume density and E* = 4/c)*/i* (van der Kruit 1988). Eliminating p*, then in terms of 
measured quantities, o*,;^ = \/2 tt G E* ft* (van der Kruit 1988) and 



Finally, we assume a flxed ratio o*^^ = 0.6 o,,,. to relate the radial and vertical velocity dispersions, which is reasonable 
for most late- type galaxies based on the limited available evidence (e.g., Shapiro et al. 2003). 

These assumptions yield disk-averaged Qstars (Equation 15) mostly in the range ~ 2-4, in reasonable agreement with 
estimates in the Milky Way (Jog & Solomon 1984; Rafikov 2001) and the expectation that stellar disks are marginally 
stable against collapse, Qstars ^ 2 (Kregel & van der Kruit 2005, and references therein). Our fixed flattening ratio 
yields nearly identical results to the fit used by Blitz & Rosolowsky (2006) to derive /i* from The scaling between 
o* and maximum rotation velocity observed by Bottema (1993) and Kregel & van der Kruit (2005) yields roughly 
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Fig. CI. — (left) Motivation for our adopted gas velocity dispersion, o"gas. We plot the median and Icr range of Cgas over the outer 
disk (0.5—1.0 r25) of each THINGS galaxy (including galaxies that are not part of this study) as a function of inclination. We see that 
variations within a galaxy are relatively small and that Cgas = 11 ± 3 km s^^ is a good description of galaxies with inclination i < 60°. 
Above this inclination, o-gas is systematically higher with more scatter as a likely consequence of projection effects in these systems. We 
adopt a fixed (Tgas (gray line) throughout this work, (right) Median if -band intensity vs. median 3.6/im intensity in 10"-wide tilted rings. 
A single scaling, I3.6 = 0.55 Ik (shown by the solid line), relates the two well. 



similar scale heights but is more sensitive to adopted structural parameters (a problem for several face-on galaxies). 
The scatter among the various methods to estimate a* or /i* and observations remains ~ 50% and this is clearly an 
area where more observations are needed (particularly measuring cr* as a function of radius, though see CiarduUo et 
al. 2004; Merrett et al. 2006). 

STELLAR SURFACE DENSITIES FROM THE IRAC 3.6 jjM BAND 

SINGS (Kcnnicutt ct al. 2003) imaged most of our sample using the IRAC instrument on Spitzer (Fazio ct al. 2004). 
Emission from old stellar photospheres accounts for most of the emission seen in the IRAC 3.6 fim band (e.g., Pahre et 
al. 2004), although we note that there may be some contribution from very hot dust and PAH features. Therefore we 
use these data to estimate radial profiles of the stellar surface density, S*. To convert from 3.6/im intensity, Is.e, to S* 
we apply an empirical conversion from 3.6^m to /C-band intensity and then adopt a standard /C-band mass-to-light 
ratio. 

We work with median profiles of /a.e, taken over a scries of 10"-widc tilted rings using the structural parameters 
in Table 4. Real azimuthal variations, e.g., due to bars or spiral arms, are lost. This is balanced by three major 
advantages from the median: 1) we avoid contamination by hot dust or PAH emission near star forming regions, a 
potential issue with the 3.6/xm band; 2) we filter out foreground stars; and 3) we increase our sensitivity by averaging 
over the ring. The first advantage avoids a serious possible bias due to confusing S* and Ssfr- The latter two allow 
us to measure T,^ out to large radii. 

To calibrate the ratio of J3.6 to if -band intensity, Ik, we compare I^e profiles to Ik profiles from the 2MASS Large 
Galaxy Atlas (LGA Jarrett et al. 2003). The profiles from the LGA are not sensitive enough to reach > r25 in most 
cases, but they yield sufficient data to measure a typical /if-to-Js 5 ratio. The right panel in Figure CI shows this 
measurement. Wc plot Ik as a function of Js.g; each point gives median intensities in one 10"-wide tilted ring in one 
galaxy. The sohd line shows a fixed ratio Js.e = 0.55 Ik (both in MJy ster~^), which matches the data very well. This 
agrees with results from Oh et al. (2008), who investigated the /r-to-3.6/im ratio using stellar population modeling 
and found only very weak variations. 

To convert Ik to we apply a fixed K-hand mass-to-light ratio, = 0.5 M0/L0 i^. This is near the mean 
expected for our sample: applying the Bell et al. (2003) relation between B — V color and mean T^, wc find Tf^ = 0.48- 
0.60 Mq/Lq^k (using global B — V colors and assuming a Kroupa 2001, IMF to match our star formation rate). This 
small range in mean motivates our decision to adopt a constant value. 

With our _ft'-to-3.6/im ratio, = 0.5 M0/L0.X, and the if-band magnitude of the Sun = 3.28 mag (Binney & 
Merrifield 1998), the conversion from 3.6 /um intensity to stellar surface density is 

S* = l-y—) cosi /s.e = 280 cosi /s.e, 

\^3.6/ 

with in Mq pc~^ assuming a Kroupa (2001) IMF and Is.e in MJy ster~^. 
The major uncertainty in Equation CI is the mass-to-light ratio, which depends on the star formation history. 



(CI) 
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metallicity, and IMF. The mass-to-hght ratio varies less in the NIR than the optical but it does vary, showing 
0.1 dex scatter for redder galaxies and 0.2 dex for bluer galaxies (Bell & de Jong 2001; Bell et al. 2003). Because 
metallicity and star formation history exert different influences on galaxy colors and T^, these variations are not 
readily inferred from colors (unlike in the optical, e.g.. Bell et al. 2003). 

In their analysis of the THINGS rotation curves, dc Blok et al. (2008) also derive from I^.e- They use J — K 
colors from the 2MASS LGA to estimate variations in T^. Their Figure 21 compares our integrated masses to those 
that they derive using color-dependent for both a (Kroupa 2001) and "diet Salpeter" (see Bell & de Jong 2001) 
IMF. Because they use the Bell & de Jong (2001) results, which have a fairly strong dependence on NIR color, they 
find Tf ~ 30-40% higher than we do in massive (red) spiral galaxies, even for matched (Kroupa 2001) IMFs. 



We combine GALEX FUV and Spitzer 24/im maps to estimate the star formation rate surface density, Ssfr, along 
each line of sight. FUV maps show mostly photospheric emission from O and B stars and thus trace unobscured 
star formation over a timescale of tfuv ~ 10-100 Myr (e.g., Kennicutt 1998b; Calzetti et al. 2005; Salim et al. 2007). 
Emission at 24/xm originates from small dust grains mainly heated by UV photons from young stars. It has been shown 
to directly relate to ongoing star formation over a timescale of T24 ~ 10 Myr (e.g., Calzetti et al. 2005; Perez-Gonzalez 
et al. 2006; Calzetti et al. 2007). We adopt this tracer because: 1) the resolution and sensitivity of the GALEX FUV 
and Spitzer 24/im maps arc both good (and well matched), 2) these data are available for our whole sample, and 3) 
the combination is directly sensitive to both exposed and embedded star formation. 

In this section, we take a practical approach, calibrating our tracer by comparing it to other estimates of Ssfr. For 
a more thorough discussion of the relationship between extinction, UV, and IR emission, we refer the reader to, e.g., 
Calzetti et al. (1995); Buat et al. (2002); BeU (2003); Cortese et al. (2006); Boissier et al. (2007). Our tracer builds 
mainly on two recent results: 1) for entire galaxies, Salim et al. (2007) showed that FUV emission can be used to 
accurately measure star formation rates (with typical rpuv ~ 20 Myr) if extinction is properly accounted for and 2) 
Calzetti et al. (2007) and Kennicutt et al. (2007) demonstrated that 24/im data could be used to accurately correct 
Ha for extinction. We combine these results using a method similar to that of Calzetti et al. (2007): via comparisons 
to other estimates of extinction-corrected Ssfr, we derive a linear combination of FUV and 24/xm intensity that we 
use to estimate Ssfr, 



unobscured SFR using the FUV-to-SFR calibration found by Salim et al. (2007); the second term measures embedded 
SFR from 24/tm and is 30% higher than the matching term in the Ha+24/im calibration of Calzetti et al. (2007). The 
additional weight reflects the fact that FUV is more heavily absorbed than Ha. More detailed motivation for Equation 
Dl is given in Appendix D.2. 

Following Calzetti et al. (2007), Equation Dl assumes the defauh initial mass ftmction (IMF) of STARBURST99 
(Leithercr et al. 1999), the broken power law given by Kroupa (2001) with a maximum mass of 120 M0. This yields 
Ssfr a factor of 1.59 lower than a 0.1-100 Mq Salpeter (1955) IMF (e.g. Kennicutt 1989, 1998a). Our FUV term is 
Equation 10 from Salim et al. (2007) divided by this value (1.59); the calibration is the same found for the Chabrier 
(2003) IMF over the range 0.1-100 Mq (their Equations 7 and 8). 



GALEX NGS FUV Maps 

We use FUV maps obtained by the GALEX satellite (Martin et al. 2005) as part of the GALEX Nearby Galaxies 
Survey (NGS, Gil de Paz et al. 2007). The GALEX FUV band covers A = 1350-1750 A with a resolution of 5.6" and a 
1.25° diameter field of view. These maps have excellent sensitivity and well-behaved backgrounds over a large field of 
view. GALEX simultaneously observes in a near- UV (NUV) band (A = 1750-2750 A). We use these data to measure 
UV colors and to identify foreground stars. 

We correct the FUV maps for Galactic extinction using the dust map of Schlegel et al. (1998). We subtract a small 
background, measured away from the galaxy. We identify and remove foreground stars using their UV color: any pixel 
with a NUV to- FUV intensity ratio > 15 (varying ±5 from galaxy to galaxy) that is also detected in the NUV map 
with > 5(7 significance is blanked. In convolution to our working resolution, blank pixels are replaced with the average 
of nearby data. We also blank a few regions with obvious artifacts. These include bright stars (e.g., in NGC 3198 and 
NGC 6946) that are usually beyond the optical radius of the galaxy and M51b, the companion of M51a. 



We use maps of 24/xm emission obtained as part of the SINGS Legacy program (Kennicutt et al. 2003) using the 
MIPS instrument on Spitzer (Rieke et al. 2004). Gordon et al. (2005) describe the reduction of these scan maps, 

which have 6" resolution. The sensitivity and background subtraction are both very good, and it is typical to find Scr 
emission at ~ r25 in a spiral. The MIPS PSF at 24/im is complex at low levels, but our working resolution of ~ 20" 
makes this only a minor concern. 
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NGC 3077, NGC 4214, and NGC 4449 arc not part of SINGS. For these galaxies we use 24/xm (and IRAC) maps 
from the Spitzer archive. We use the post basic calibrated data produced by the automated Spitzer pipeline. 

As with the FUV maps, we subtract a small background from the 24 /um maps, which we measure away from the 
galaxy. We blank the same set of foreground stars as in the FUV maps. In convolution to our working resolution, 
these pixels are replaced with the average of nearby data. We also blank the edges of the 24/xm maps perpendicular 
to the scan, which are noisy (and outside the optical radius) and the same artifacts blanked in the FUV maps. 

SINGS Ha 

The SINGS fourth data release includes Ha maps, which we use to compare Ssfr derived from Ha, FUV, and 
24/im emission in 13 galaxies. We convert Ha to Ssfr following Calzetti et al. (2007) and the SINGS data release 
documentation. We correct for [Nil] contamination following Calzetti et al. (2005) and Lee (2006) and correct for 
Galactic extinction using the Schlegel et al. (1998) dust maps. We check the flat-fielding by eye and mask regions or 
flt backgrounds where necessary. 

Motivation for the FUV+24 to EgpR Relation 

Because 24/im lies well short of the FIR peak for a typical galaxy SED (e.g., Dale & Helou 2002), a measurement in 
this band does not directly trace the total IR luminosity. Therefore using 24^m to measure embedded SFR relies on 

modeling of the IR SED or empirical calibration against other estimates of extinction. Calzetti et al. (2007) compared 
Ha and 24/Ltm to Paschen a (Paa) emission, a tracer of ionizing photons largely unaffected by extinction. They showed 
that a linear combination of Ha and 24/im, 

SFRxot = SFR^"„°'^^™''^'* (Ha) + SFR^;^'^^'*^^'* (24^m) 

SFRTot = 5.3 X lO"'^^ (i [Ha) + 0.031 L (24/xm)) , (D2) 

matches SFRxot inferred from Paa for 220 individual star forming regions in 33 nearby galaxies and that the same 
calibration also works well when integrated over a large fraction of a galaxy disk. Here L (Ha) is the luminosity of Ha 
emission from the region in erg s~^ and L (24/xm) = i'24^in (24/xm), also in erg s~^, is the specific luminosity of the 
region times the frequency at 24/im. SFR is the star formation rate in that region in Mq yr~^. 
We require an analogous formula to combine FUV and 24/im data, 

SFRTot = SFR^'Jj°^«^"'-<'<i (FUV) + SFR^'S^^'^^'^'^ (FUV, 24/xm) . (D3) 

The first term is the SFR implied by a particular FUV luminosity taking no account of internal extinction. The second 
term is the SFR that can be attributed to FUV light that does not reach us — i.e., the extinction correction for the 
first term — , which we infer from the 24/im luminosity and may also depend on the ratio of FUV to 24/im intensity. 

We adopt the first term in Equation D3 from Salim et al. (2007), who studied the relationship between FUV emission 
and SFR in ~ 50, 000 galaxies, combining multi-band photometry with population synthesis modeling and comparing 
to Ha emission. They found 

SFR^'Jj"^^'^"''^^ = 0.68 X 10-2^ (FUV) , (D4) 

with SFR in Mq yr~^ and L^, (FUV) in erg s~^ Hz~^. This yields SFRs ~ 30% lower than the relation given by 
Kennicutt (1998b) because of metallicity, model, and star formation history differences between their sample and 
Kennicutt's model. 

We calibrate the second term in Equation D3 in two ways: 1) we use simple assumptions to extrapolate 
SFRI'J^'^^d'i'^d (FUV,24/im) from SFRfi'^'''''^'*'''^ (24/xm), which was measured by Calzetti et al. (2007); and 2) we make 
several independent estimates of SFRxot — by comparing SFRxot with FUV emission, we directly measure the second 
term in Equation D3. We phrase both analyses in terms of the factor Wpw, defined as 

qp'D embedded 

The numerator is the second term in Equation D3 and the denominator is the relation between embedded Ha and 
24/im emission measured by Calzetti et al. (2007, Equation D2). To measure Wpuvi we combine Equations D3 and 
D5 to obtain 

SFRxot = SFR^'ff^^^"''^'^ (FUV) + Wfuv (FUV, 24/im) SFRg^'"''^^'''^ {24fxm) (D6) 
and solve for W^fuv in terms of measurable quantities 

^ ^ SFRTot - SFRru°r--'-^FUV) 
"""^ SFRg^^^'^'i"'* (24/xm) ' ^ ' 

So that to estimate Wpuv over a line of sight we require FUV and 24/im intensities and an estimate of SFRxot- 
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Simple Extrapolation 

In conjunction with a direct measurement, it is helpful to have a basic expectation for VFpuv- We calculate this by 
combining a Galactic extinction law and a typical nebular-to-stellar extinction ratio. In terms of Ha extinction, ^Hai 
and FUV extinction, ApuV) Equations D2 and D3 are 



SFRxot = SFRS"^"''''^"'"^'* io^Ho/2.5 
SFRxot = SFR^^°^'^^"'''"^ io^Fuv/2.5 ^ 

For a Galactic extinction law ^puvMii ~ 8.24/2.33 (Cardelli et al. 1989; Wyder et al. 2007). We may also expect that 
FUV originates from a slightly older and more dispersed population than Ha. If we assume a typical nebular-to-stellar 

extinction ratio of Ahq/^a ~ 2 (Calzctti ct al. 1994; Roussel et al. 2005), then wc expect Apuv/^Ha ~ 1-8 (if FUV 
comes mostly from a very young population coincident with Ha, we instead expect Apuv/^Ha ~ 3.6). Combined with 
Equations D3 - D8 these assumptions yield 



+ 1 = f ^ + l) where / = SFRj^-^r j; (FUV) 

Z + TUfuv-I V ./ / SFRg^'^^'i'i^d (24Aim) ' ^ ' 

which wc may solve for Wpuv given /, the ratio of observed FUV-to-24y^m intensities (in SFR units). 

For Ana = 1.1 mag, a typical value in disk galaxies (Kcnnicutt 1998b), / ~ 0.26 and Equation D9 suggests 
Wpuv ~ 1-3. For higher ^Ha, expected for the inner parts of spiral galaxies or arms, / will be lower and we expect 
lower values of PFpuvj approaching Wpuv = 1 where both FUV and Ha are almost totally absorbed (and SFRxot is 
determined totally from 24/im emission). For lower Ahq, e.g.. expected in dwarf galaxies or the outer parts of spirals, 
we expect Wpuv to approach the ratio of extinctions, 1.8, in the optically thin case. 

Measuring TUpuv 

We measure VFpuv directly from observations by comparing FUV and 24/im emission to various estimates of SFRxot- 
We perform these tests in the 13 galaxies with SINGS Ha data. Over a common set of lines of sight where we estimate 
Ha, FUV, and 24/im to all be complete, we estimate Esfr and Wpuv (from Equation D7) in 5 ways: 

1. Combining Ha-|-24/xm using Equation D2 (Calzetti et al. 2007). 

2. Prom 24/im emission alone, using the (nonlinear) relation found by Calzetti et al. (2007, their Equation 8). 

3. From Ha alone, taking A-na = 1-1 mag, a typical extinction averaged over disk galax;ies, though not necessarily 
a good approximation for each line of sight (Kennicutt 1998b). 

4. From Ha emission, estimating Ana from Shi and Sh2 following Wong & Blitz (2002). We assume a Galactic 
dust-to-gas ratio and treat dust associated with H I as a foreground screen obscuring the Ha while treating dust 
associated with H2 as evenly mixed with Ha emission. 

5. From FUV emission, estimating Apuv for every line of sight applying the relationship between FUV-to-NUV 
color and ^puv measured for nearby galaxies by Boissier et al. (2007). 

In principal, the first method is superior to the others because Calzetti et al. (2007) directly calibrated it against 
Paa, and because it incorporates both Ha and IR emission, offering direct tracers of both ionizing photons and dust- 
absorbed UV light. The other four methods offer checks on SFRxot that are variously independent of 24/im, FUV, or 
Ha emission, allowing us to estimate the plausible range of both W^puv and the uncertainty in Sspr. 



Derived Relation 

Figure Dl shows the results of these calculations. In the top panel, we plot the normalized distribution of Wpuv for 

each estimate of SgpR . The bottom left panel shows how each distribution of Wpuv depends on the FUV-to-24/im 
ratio, / (Equation D9). The bottom right panel shows how lUpuv varies with normalized galactocentric radius. 

The median T4^fuv derived in various ways spans a range from ^1.0 1.8. The two 24/im based methods (blue 
and gray) both yield Wpuv ~ 1.3 with relatively narrow distributions. Using FUV and UV-colors yields the highest 
expected Wpuv; ^ 1-8; estimating Ayi^ from gas yields the lowest WpuV) peaked near ~ 1.0, though the distribution 
is very wide. This range of values agrees reasonably with our extrapolation (seen as a dash dotted curve in the top 
right panel), which also lead us to expect a typical T4^fuv of 1.3 and a reasonable range of 1.0 1.8. 

The bottom panels show that while individual methods to estimate M^fuv do exhibit significant systematics (partic- 
ularly at very high and low of /), simply fixing I¥puv = 1-3 is a reasonable description of most data (the dashed lines 
in the center panel bracket 80% of the measured /). Wpuv does not have to be constant. Indeed, we expect it to 
vary with / based on simple assumptions and very basic arguments. However, a constant Wpuv is consistent with the 
data and is also the simplest, most conservative approach. Therefore, this is how we proceed: taking VFpuv = l-^^oJ, 
Equation D3 becomes 
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Fig. D1. — Wpuvi the calibration of the term to estimate the SFR from a Unear combination of FUV and 24/im emission. We 

measure Wpuv pixel-by— pixel by comparing FUV and 24/im intensity to EsFR estimated in five ways: (blue) combining Ho and 24/im; 
(gray) using only 24/.tm; (red) using Ho, estimating extinction from the gas; (green) using Ho, assuming a typical extinction; and (purple) 
using FUV emission, estimating Apyy from the UV color. We plot the resulting Wpuv ™ three ways: (top) as normalized histograms; 
(bottom left) as a function of /, the ratio FUV to 24/im emission along a line of sight (see Equation D9); and (bottom right) as a function 
of galactocentric radius normalized by r25. The hatched regions in the bottom panels show the median trend ztlo" for each case. In each 
panel, we indicate our adopted IVpuv = 1-3^q'|. The dash— dotted curve in the bottom left panel shows the expectation for a typical 
extinction law and nebular-to-stellar extinction ratio and the vertical dashed lines show the range of / that includes 80% of the data. 





Fig. D2. — Five estimates of SgpR (j/-axis) as a function of SgpR predicted by our combination of FUV and 24fii[a. The color scheme is 
the same as Figure Dl and the methodology used to derive Ssfr for comparison is labeled in each plot. In the bottom right panel, each 
point shows the integrated SFR for a galaxy derived from Ha as a function of the SFR derived from FUV+24^m emission. In all plots, 
solid lines show slopes of 0.5, 1, and 2. 



SFRTot = 0.68 X IQ-^^ (FUV) + 2.Utoil x 10"-*^ L (24^ni) , (DIO) 

We convert Equation DIO from luminosity to intensity units using J^24/^m = 1-25 x 10^'^ Hz, 1 MJy= 
10-1"^ erg Hz-i cm ^, and L,y = Air A I^, where A is the physical area subtended by the patch of sky being 
considered. This yields Equation Dl, 



EsFR = 8.1 X 10-2/fuv + x 10-^/24 , (Dll) 



with /puv and in units of MJy ster ^ and Ssfr in units of Mq kpc yr 



-0.7 

2 ,,^-1 



Uncertainty in SgpR 

In Figure D2, we plot our 5 alternate estimates of Ssfr (y^axis) as a function of SgpR derived from Equation Dll 
(x-axis). Each point corresponds to a 10"-wide tilted ring. In the bottom right panel, we plot the SFR integrated 
over the disk (over rgai < ^25) as a function of SFR estimated from nebular line emission by Kennicutt et al. (2003) 
and Ha by Lee (2006). Solid lines in all six panels show the line of equality plus or minus a factor of 2. 

If we adopt the naive tack of treating all approaches as equal, the aggregate data in Figure D2 yield a median ratio 
Ssfr (other)/SsFR(FUV + 24) « 1.05 with « 0.22 dex (i.e., ~ 65%) la scatter. The dominant sources of this scatter 
are the choice of "other" Ssfr and galaxy-to-galaxy variations. Once a galaxy and methodology are chosen, the data 
tend to follow a fairly well-defined and often nearly linear relation. For comparison, the 24/im part of the Ha+24/im 
calibration has « 20-30% uncertainty (Calzetti et al. 2007; Kennicutt et al. 2007) only considering star forming peaks. 
In light of the wider range of star formation histories and geometries encountered working pixel-by-pixel or averaging 
over whole rings, the estimate of ~ 65% seems quite reasonable. Comparing our integrated SFRs (Figure D2, bottom 
right) to those estimated by IIHUGS Lee (2006) and SINGS Kennicutt et al. (2003) bears out this estimate; we match 
these estimates with a similar scatter. Another view of this comparison may be seen in Appendix F, where we present 
radial profiles of Ssfr based on Ha on the same plots as our FUV+24/im profiles. 

Despite the overall good agreement between our Ssfr and other estimates. Figures Dl and D2 do show systematic 
differences among tracers. We note several of these before moving on: 
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1. Using only 24/im emission (gray) yields a low estimate of SgpR for the two low mctallicity galaxies in our 
comparison sample: Holmbcrg II and IC 2574. Dust is known to be deficient in these galaxies (Walter et al. 
2007), which is likely to lead to a breakdown in the fit between 24/xm emission and Paa. This effect, already 
recognized by Calzetti et al. (2007), highlights the importance of including a non-IR component in a SFR tracer. 

2. Estimating Aua from gas (red, Wong & Blitz 2002) yields very high Ssfr (and high Wpuv) in the inner parts 
of galaxies. This underscores the complexity of the geometry and timescale effects at play; it is extremely 
challenging to reverse engineer the true luminosity of a heavily obscured source knowing only the amount of 

nearby interstellar matter. These high values are almost certainly overestimates; stellar feedback, turbulence, or 
simply favorable geometry likely always allows at least some light from deeply embedded H II regions to escape. 

3. Particularly at low Ssfr, inferring Apuv from UV colors (purple) yields higher embedded SFR than using 24/zm 
emission (and this method appears to completely fail in NGC 6946, the horizontal row of points). A possible 
explanation is that where Ssfr is relatively low, the UV originates from a somewhat older (and thus redder) 
population (e.g. Calzetti et al. 2005); the FUV-UV color relation depends on the recent star formation history 
(e.g., differing between starbursts and more quiescent galaxies, Boissier et al. 2007; Salim et al. 2007). 

This discrepancy (and the close association between our SFR tracer and stellar mass seen in the main text) 
argues for a comparison among mctallicity, stellar populations, and mid-IR emission that is beyond the scope 
of this paper. We restrict ourselves to a first-order check: we compare the ratio of 24/im to 3.6/im and FUV 
to 3.6/im emission in our sample to those in elliptical galaxies, which should be good indicators of how much 
an old population contributes to 24/xm or FUV emission. Very approximately, in ellipticals I2A/H.Q ~ 0.1 (Temi 
et al. 2005; Dale et al. 2007; Johnson et al. 2007), with ~ 0.03 expected from stellar emission alone (Helou 
et al. 2004), while /fuv/^3.6 ~ 2-4 x 10~^ (Dale et al. 2007; Johnson et al. 2007, taking the oldest bin from 
the latter). We measure 124/13.6 Ipw/Iss for each ring in our sample galaxies and compare these to the 
elliptical colors. In both cases, only ~ 5% of individual tilted rings have ratios lower than those seen in elliptical 
galaxies and the mean color is ~ 10 times that found in elliptical galaxies, though the ratio /fuv/^3.6 shows 
large scatter due to the effects of extinction. Both the 24/im and FUV bands do appear to be dominated by a 
young stellar population almost everywhere in our sample. Discrepancies among various tracers thus seem likely 
to arise from the different geometries and age sensitivities of FUV (r ~ 100 Myr), Ha (r ^ 10 Myr), and 24/im 
(likely intermediate) emission. 

Finally, we emphasize that uncertainties inferred via these comparisons mainly reflect the ability to accurately infer 
the total UV light or ionizing photon production from young stars. They do not include uncertainty in the IMF, 
ionizing photon production rate (e.g., at low metallicity) , or any of the other factors involved in converting an ionizing 
photon count or FUV intensity into a SFR. 

RADIAL PROFILES 



Table El presents radial profiles of Sjji, 5]h2, S*, and Ssfr- Combined with kinematics, which may be calculated 
using Equation Bl taking Vdat and Zflat from Table 4, these profiles are intended to provide a database that can be 
used to test theories of galaxy-wide star formation or to explore the effects of varying our assumptions. Results for all 
galaxies are available in an electronic table online. Table El in the print edition shows the results for our lowest-mass 
spiral galaxy, NGC 628, as an example. 

The individual columns are as follows. Ring identifiers: (1) galaxy name; galactocentric radius of ring center (2) 
in kpc and (3) normalized by r25- Mass surface densities (in Mq pc~^) along with associated uncertainty of (4-5) 
Hi; (6-7) H2; and (8-9) stars. Star formation rate surface density, Ssfr, with associated uncertainty (10-11) from 
combining FUV and 24/zm emission in units of 10~^ Mq yr~^ kpc~^; and the individual contributions to Ssfr from 
(12) FUV and (13) 24/im emission (i.e., the left and right terms in Dl) in the same units. 

We derive radial profiles from maps using the mean (for Shi, 5]h2, SgpR.) or median (S*) value within 10"-wide 
tilted rings (so that the rings are spaced by half of our typical working resolution) . The rings use the position angle 
and inclination in Table 4, adopted from Walter et al. (2008). We adopt the THINGS center for each galaxy (Walter 
et al. 2008; Trachternach et al. 2008) except for Holmberg I, where we use the dynamical center derived by Ott et 
al. (2001) rather than the photometric center. We consider only data within 60° of the major axis, measured in the 
plane of the galaxy. This minimizes our sensitivity to the adopted structural parameters, which most strongly affect 
the deprojection along the minor axis. Where there are no data, we take Ehi = and Sh2 = 0. These are regions 
that have been observed but masked out because no signal was identified. We ignore pixels with no measurement of 
5^SFR, these are simply missing data. 

We take the uncertainty in a quantity averaged over a tilted ring to be 

^^^^ (El) 

Y -^piXjringZ-^piXjbeam 

where ctrms is the RMS scatter within the tilted ring, A^pix.ring is the number of pixels in the ring, and iVpix,beam is 
the number of pixels per resolution element. This a captures both random scatter in the data and variations due to 
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^The full table of radial profile data is available in electronic form. 



azimuthal structure within the ring. It does not capture systematic uncertainties, e.g., due to choice of Xco or star 
formation tracer, discussed in these appendices. 

ATLAS OF MAPS AND PROFILE PLOTS 

In Figure F, we present maps, profiles, and calculations for individual galaxies. Each page shows results for one 
galaxy. The top row shows maps of atomic gas (Shi), molecular gas (Eh2), and total gas (Sgas = Shi + Sh2)- The 
second row shows unobscured (FUV), dust-cmbcdclcd (24/im), and total star formation surface density (SgpR,). These 
maps use a color scheme based on the modified magnitude system described by Lupton et al. (1999); a bar to the 
right of each row of plots illustrates the scheme. The gas maps and star formation maps for each galaxy use a single 
color scheme, but the scheme does vary from galaxy to galaxy, so care should be taken when comparing different 
galaxies. Also note that we construct the table to show empty values below our working sensitivity (i.e., any data 
below Sgas = 1 M0 pc~^ or Esfr = 10~^ Mq yr~^ kpc~^ appear as white) but the data (especially THINGS) often 
show evidence of real emission below this value. We refer the readers to the original data papers for more information 
on each data set. 

The dotted circle indicates the optical radius, r25, in the plane of the galaxy for the structural parameters given in 
Table 4. A small black circle in the bottom right panel shows our working resolution. 

In the left panel on the third row, we plot mass surface density profiles. We show H I (blue), H2 (magenta, where 
available), stars (red stars), and total gas (thick gray profile). Vertical dotted lines indicate 0.25, 0.5, 0.75, and 1.0 
times r25. Horizontal dotted lines show fixed mass surface density. The thick gray vertical line shows where the 
intensity scale for the images is set, i.e., 0.1 r25. 

In the right panel on the third row, we plot star formation rate surface density profiles. We show the total Ssfr 
(thick gray profile) and the separate contributions from dust-embedded (green, 24/im) and unobscured (blue, FUV) 
star formation, which add up to Ssfr- Where they are available, we plot Ssfr from the SINGS DR4 Ha (red) and 
points measured from the Ha profiles of Martin & Kennicutt (2001) (magenta) and Wong & Blitz (2002) (purple). 
All Ha profiles assume 1.1 mag of extinction (a typical average value in disk galaxies, Kennicutt 1998b). The other 
markings are as in the left panel. 

In the left panel of the fourth row, we show the observed SEE for the galaxy. We use the same color scheme as in §4, 
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i.e., magenta points indicate rings where Sh2 > Shi, blue points show rings where Sh2 < Shi, and red arrows indicate 
upper limits. The ensemble of points in this panel combine to form Figure 1. Dashed, dotted, and dash-dotted lines 
show the SFE predicted following the method described in §6 with no threshold applied (the thresholds appear in the 
right panel). The other markings are as in the panels on the third row. 

In the right panel of the fourth row we show azimuthally averaged values for thresholds described in §2.2. We expect 
widespread star formation (conditions are "supercritical" ) where the value of a profile is below 1 (the shaded area) 
and isolated or nonexistent star formation (conditions are "subcritical" ) above 1. We plot the Toomre Q parameter 
for a gas disk, Qgas (black), and for a gas disk in the presence of stars, Qstars+gas (green). We show the shear criterion 
described by Hunter et al. (1998a), Scrit.A/Sgas in purple and the condition for the formation of a cold phase given by 
Schaye (2004), Sso4/Sgas in orange. The other markings are as in the panels on the third row. 
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Fig. F1. — Atlas of data and calculations for DDO 154. 
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Fig. F. — Atlas of data and calculations for NGC 4214. 
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Fig. F. — Atlas of data and calculations for NGC 2976. 
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Fig. F. — Atlas of data and calculations for NGC 3077. 
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Fig. F. — Atlas of data and calculations for NGC 2403. 
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Fig. F. — Atlas of data and calculations for NGC 628. 



58 



Leroy et al. 



NGC 3198 



u 
m 

Ll 

t- 

nS 

Ha 



m 
n 
(J 
I- 

•a 
< 



300 


300 
D 

-soo 



' 

V \ 

■ill' 

1 t ■ 1 ■ ' ' 1 1 ■ t 1 ■ . t 


\ ♦ 

*, • * 
\ '* \ 

'a \ 

'hi ■ 

■ 1 1 1 ' 1 ' 1 ■ ■ 1 1 1 1 ■ 


Tolat Caa 
' ll • 

— , . 1 . . . 1 — , 






Embedded SPH 


Total sm 
















% 














'Hi » 






■ 


K 

, 1 , ... 


m 

I (.11 






(arcae-c) 


2Q0 -EOP 
AoL (arL'acL"] 



ft. 



]0 



10' 



■1 

r — I 

T 




Fig. F. — Atlas of data and calculations for NGC 3198. 
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NGC 3184 
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Fig. F. — Atlas of data and calculations for NGC 3184. 
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NCC 4736 
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Fig. F. — Atlas of data and calculations for NGC 4736. 
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NGG 3351 
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Fig. F. — Atlas of data and calculations for NGC 3351. 
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NGC 6946 
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Fig. F. — Atlas of data and calculations for NGC 6946. 
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Fig. F. — Atlas of data and calculations for NGC 3627. 
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NGC 5194 
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Fig. F. — Atlas of data and calculations for NGC 5194. 
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Fig. F. — Atlas of data and calculations for NGC 3521. 
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NGC 2941 





Fig. F. — Atlas of data and calculations for NGC 2841. 
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Fig. F. — Atlas of data and calculations for NGC 5055. 
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NGC 7331 
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Fig. F. — Atlas of data and calculations for NGC 7331. 



